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Any  one  can  be  expressed  in  terms  of  two  others 
Deformations  can  be  decomposed  into  two  types 
Hydrostatic  compression:  constant  shape 
Shear:  constant  volume 
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Wave  Equation  -  No  Dispersion 
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From  Conservation  of  Mass  (7a) 


Qa  _  _  i-  9p<? 


Conservation  of  Momentum 


Q©  - 


r° 


First  Law  -  Conservation  of  Energy 


C.  o (SUb 
o 


Assume 

-  T><0  ^ 


■  tart.  --W) 


Solution  gives 

Cc^.-  fl) 


CiO 


Pure  Real 

Classical  Absorption 


24 


Wave  Equation  -  No  Dispersion 


From  Conservation  of  Mass  (7a) 
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Classical  Absorption 


Eq.  2  Eq.  12 

FigTl 


Assume  p  f  3 


2ju 
S>x  * 


which  gives 

a--Uv  e-^  e  Wa) 


where 

£  a)/C0  oL  =  £lj  (  Cq  ) 
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(_22.) 


'  For  water,  <=*  ^  lV  ^ 


Classical  Absorption  Continued 
Heat  Conduction 


Fourier’s  Law 


Some  approximations  give 

~  l cs'-o 

dp  po  Cq 


-  j 

Note  reappearance  of  co  term  ; 
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Relaxation  Absorption/Dispersion  Modes 


Energy  Modes 

Translation,  Rotation,  Vibration 
Collision  Frequency  ~10n/second  in  gas  at  STP 
Translational  Relaxation  7^-  w?  7C 
Rotational  Relaxation 


ini  CHO  ''Sr0 
-iM  (Ml)  v  3  cr  ‘i 
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Semi-Classical 

Trajectories 
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b.  _ 


'  X 


Due  to  energy  spacing 

i 


29 


\f-r 
v-i R 

v  —\f 
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Summary  of  Absorption  Mechanisms 

Viscosity  Momentum 

Thermal  Conduction  Energy 


Internal  Relaxation 


Internal  Energy 
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Relaxation  Absorption 


The  Fun  Part! 


Experimental  Techniques 


Experimental  Considerations 
Viscous  and  Thermal  Losses  in  Gas  a  f2 

Viscous  and  Thermal  Losses  at  WaIIsaf1/2 

Radiation  or  Leakage  -  Greater  at  Low  Frequencies 

Spreading  Losses  -  6dB/doubling 

Rotational  Relaxation  a  f2 
=  Tough  to  measure  especially  in  liquids 


Figure  3.  Schematic  diagram  of  the  apparatus  for  measuring  the  attenuation  of  sound. 


Fig.  4.  Diagram  of  the  experimental  system  for  measuring  sound  absorption. 


Spectrophone 


Pulsed  Laser  Test  Cell 


Figure  2  shows  typical  curves  for  absorption  per  wavelength  and  velocity  dispersion  to  a  single 
relaxation  process.  The  example  here  is  Fh  at  102°C.  The  figure  compares  curves  representing  the  above 
theory  with  measured  values.  ^ 
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Fig.  5.  Relaxation  absorption  coefficient  per  wavelength  vs  log  (frcquency/pressure)  for  chlorine.  The 
solid  curves  are  theoretical  absorption  with  values  of  Am  and  fm  adjusted  to  give  the  best  fit  of  the 
experimental  points. 


40 


Conclusions  from  Experiments 

log  P10  a  T 1/3 

P10  very  sensitive  to  AE 

P10  very  sensitive  to  At 

P10  depends  upon  geometry  of  the  mode 

P10  depends  upon  interaction  potential 
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Relaxation  Studies  in  S02 

£*«rjy  U»«l 


50t  Oi 


v 


Conclusions 

Argon  1/10  &  02 1/4  as  effective  as  S02  bend 


v-v  Coupling  with  1361  cm'1  between  S02  and  02 


v-t  Transfer  increases  most  rapidly  with 
temperature 
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Relaxation  in  CS2 

Gas  fmax  -370  ■  560  KHz/Atm 

Liquid  fc  -1013/sec 

fmax ~  10s/  sec 
Tj  =  22ns 
t2  =  30ns 
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Absorption  of  Sound  in  Air 
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Fig.  8.  Components  and  general  behavior  of  total  air  absorption  in  air  in  terms  of  (a)  loss  per 


wavelength  and  (b)  loss  per  unit  distance: 


la)  /i  ■ 


:t  i.o 
i  +  ( fO: ' 


(b)  o  - 


(f.01 
I  +  </.£r 
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DIFFUSION 

ABSORPTION 


“ON  THE  WAVING  OF  HANDS” 
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Diffusion 


*  4T>Nn 


L  wT 


QfrOjr .  -i* 


Let  M2  -  Mt  get  large 
Dusty  Gas 
Biot  Solid 


Boltzmann’s  Transport  Eq. 

Transport  described  by  distribution  function 
f(c,r,t)  and  its  time  derivative 


Collsjtch-* 


Collisions 


Collisions  take  f  from  f  to  f 


+  C3 


t-p'-p,'-  r-n ) cl«*c dc ' 


Where  i  signifies  collision  partners 

is  the  reversible  probability 
3  a  constant  f  physics  enters 

here 


Solutions  are  numerous  but  not  exact 


Fig.  10  Sound  propagation  in  n 
author  and  Academic  Press,  Inc, 


Free  Molecule  Propagation 


Assume  vibrating  surface 

/y 

00-^0  xu**  <oxt 


-f+Cx=t>)  $  -f_  ^=0^ 


Cwl 


mean  velocity  I 

of  molecules  absorption 

proportional 
to  distance 


phase  velocity 
increases  with 
distance 


At  large  distances, 

only  fast  molecules  support  propagation 
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TIME  (ms) 


Fig.  1 1 .  Filtered  microphone  response.  In  moving  from  the  central  graph  (e),  only  one  of  the  four 
variables  changes.  For  (c),  (e),  and  (0  the  sound  frequency  increases;  (d)  is  the  fundamental  of  the  60-cm 
tube;  (e)  is  the  fundamental  for  the  30-cm  tube.  Corresponding  frequencies  are  approximately  340,  680, 
and  1 360  Hz.  For  (b),  (e),  and  (h)  the  pressure  increases  from  20  to  40  to  80  Torr.  For  (a),  (e),  and  (i)  and 
H2  concentration  increases  from  0  to  10%  to  20%.  For  (g),  (e),  and  (c)  the  energy  per  mole  deposited  by 
the  discharge  in  the  gas  increases  from  3700  to  7000  to  9000  J/mol.  If  all  of  this  energy  were  to  go  into 
vibrations,  the  corresponding  initial  vibrational  temperatures  would  be  1600,  2100,  and  2400  K. 


Energy 


Conservation 


dQ 

at 


q  =  pscudoviscous  term 

dQ  =  rate  at  which  energy  is  added  to  or  taken  from  dement 
3t 

Conductivity  of  plasma 

a  =  e2  \  pm^T) 


Xc  =  electron  mean  free  path 

4.173  x  IQ'10  (A,  +  A2)  T  IC 

2x10  (1  *  Aj)  +  A]  (ajav 

At  =  fraction  of  atoms  which  have  been  ionized 
A2  =  fraction  of  atoms  which  have  lost  second  electron 
|aiJav  =  average  electron-ion  cross  section 

=  2.8  x  10->6  T-2  [(Aj  +  A2)  /  (Ai  +  A2)]2 
X  log  {1.727  x  10“5  [(A1+A2)  /  (Aj  +3A2)]  T  (A1p)“1-/2 


Fig.  12.  A  few  equations. 
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Input 

1=  IoOT^-e-*) 


Io  =  4x  104  A 
a  =  4  x  104  see-1 
p  *4  x  105  see-1 


where  j  is  a  zone  outward  from  center  of  discharge 


Losses 


Thermal  conduction 
Radiative  Energy  Transport 

Bremsstrahlung  -  Increases  as  T 1/2 
Black-body  type  radiation  -  Increases  as  T4 


Fig.  13.  Energy  balance. 
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ENERGY  (mtV) 


FIG.  4.  spectra  of  transverse  TAj  phonons  measured  along 
(1 10).  The  polarisation  vector  is  along  IlIOl.  i.e.,  the  scattering 
plapc  is  OihQ)  (sample  MIT-1). 


soft  but  sharpens  significantly.  We  confirmed  by  means  of 
additional  measurements  in  other  Brillouin  zone*  that  thn 
mode  has  Z*  symmetry  (in  Weber  s  notation).  . 

Although  we  expect  that  the  n-t  degeneracy  is  lifted 
below  Te,  we  did  not  observe  the  expected  splitting  of  the 
modes,  either  because  the  energy  resolution  of  the  spee> 
trometer  was  too  coarse  or  one  of  the  mode*  jenormahzei 
only  weakly  and  remains  buried  under  the  superlattk* 

peak.  Further  neutron  scattering  studies  are  necessary  to 
elucidate  this  point  further.  .  . 

After  identifying  the  lattice  dynamics,  which  leads  to 
the  tetragonal-orthorhombic  phase  transition,  we  per¬ 
formed  a  detailed  study  of  the  low-lying  phonon  branches 
in  order  to  fipd  evidence  for  the  predicted  breathing-mode 
instability4  of  the  LA  mode  near  the  zone  boundary.  The 
LA  phonons  were  measured  in  crystals  MIT-1  and  MiT-2 
using  very  fine,  steps  in  momentum  in  order  to  find  the  in¬ 
stability  (see  Fig.  3  for  some  representative  scans).  The 
phonons  remain  sharp  and  the  dispersion  has  a  conven¬ 
tional  form,  Le.,  the  LA  phonon  energy  increases  mono- 
tonically  with  increasing  q.  This  result  is  in  contrast  toan 
observation  in  isostructural  LajNiOa  where  the  phonon 
energy  decreases  near  the  zone  boundary.  Moreover, 


FIG  5  (a)  Summary  at  the  low-lying  phonon  branch**  in  Lai-sSrxCuO*  (aamplea  MIT-1  and  MIT- 2).  The  hsm  are 
the  eye.  Tks  modeshave  been  labeled  according  to  Ref.  4.  The  dupemoncorve.  are  only  wttkiy 
oceplioe  of  the  TO  phooon  near  tbe  X  point  bee  inaet  (c)l.  Cb)  Calculated  dupenmn  curm.  Filled 

Meeento  (bar  phMoas)  end  the  open  circle*  indicate  the  phonon*  with  I,  *ymmetry,  which  are  renormalized  by  interaction,  wd» 
the  ronrtnrtirm  electron*.  The  phonons  with  Lt  symmetry  do  not  reoonnalize  (after  Weber  (Ref.  4)1. 


imtpmcitv 
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RUS  has.... 

Made  significant  advances  in  our  understand  of  the 
structural  phase  transition  in  High  Tc  materials 

Advanced  our  understanding  of  the  electronic  states  in 
new  narrow-gap  insulators 

Provided  extensive  data  on  the  elastic  constants  of 
many  new  materials  for  structural  and  high 
temperature  applications 

Generated  a  new  area  of  measurement  science  with  its 
own  annual  conference 

Hosted  PhD  thesis  students  and  postdocs 

Generated  many  patents  on  the  use  of  resonances  for 
non-destructive  testing 

Assisted  in  LANL  programmatic  efforts  involving 
intelligence  and  stockpile  stewardship 

Transferred  technology  to  a  private  company  that  now 
manufactures  RUS  NDT  systems  in  use  on,  for  example, 
GM  assembly  lines. 


ELASTIC  UjAtV&S 
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Why  study  ultrasound? 


AG  =  Ghot  -  Gcold 

AG  is  the  Gibbs  free  energy 
change  across  a  phase  boundary 


AV  =  4-A  G 
cP 

AV  is  the  volume  change  across 
the  boundary,  and  is  zero  for  second 
order  phase  transitions 


1 
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^  hot  ^  cold 


i  dAV  ^  i  (Fag 
v  sp  ~  V  cP2 


A*ris  the  discontinuity  in  bulk  modulus 


across  the  boundary,  which  ties  directly 


to  the  speed  of  sound,  and  which  is  not 
zero  for  a  second  order  transition 
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Signal/noise  comaparison  of  pulsed  and  resonant  measurements 


parameter 

Impulse 

Swept  Sine 

drive  power  per  unit 
bandwidth 

peak  power/full 
bandwidth 
106/109=.001 

peak  power/sweep 
rate 

1/100=0.01 

noise  bandwidth  for 
complete 

measurement  using 
optimum  receiver 

109 

number  of  modes  x 
width  of  each  mode  x 
10=104Hz 

drive  duty  cycle 
(typical) 

10'3 

1 

detect  duty  cycle 

1 

1 

square  root  of  all 
factors,  which  is  a 
measure  of  S/N 

r^ 

O 

X 

CO 

10'3 

Table  I.  Signal-to-noise  comparison  between  impulse  (pulse-echo)  and 
swept-sine  (RUS)  resonance  measurement  methods  for  a  measurement  of 
a  1  cm  sample  with  resonances  having  a  Q=104,  using  10  modes  over  0.5 
MHz-1 .5MHz=106  Hz  bandwidth  to  obtain  an  elastic  modulus.  Note  that  the 
pulse-echo  measurement  provides  about  0.1%  absolute  accuracy  at  best, 
compared  with  about  0.01%  for  the  best  RUS  measurements. 
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Instrumentation  Block  Diagram 

Transducers 

Transmitting  Receiving 


Low  temperature  RUS  cell 
with  a  2mm  rectangular 
parallelipiped  sample  set 
between  2  diamond/LiNb03 
transducers  using  no 
coupling  fluids.  The  force  on 
the  sample  is  equivalent  to  1 
gm  weight.  This  cell 
operates  between  IK  and 
400K. 


Frequency  (MHz) 


Ultimate  accuracy  determined  by  geometry-for  a  Si3N4  ball 
bearing,  geometry  errors  are  less  that  1  part  in  105.  So  are  the 
modulus  errors! 


Table  1 

Resonant  ultrasound  measurement  of  a  0.63500  cm  diameter  Si,N4  ceramic  sphere  with  a  density  of  3.2325  g/cm  .  fm  are  measured 
frequencies,  fT  are  fitted,  n  is  the  mode  number,  k  is  our  designator  (to  be  discussed  below)  for  the  symmetry  of  the  mode  and  i  is 
in  essence  the  harmonic  number  of  each  symmetry  type.  Multiple  entries  indicate  the  mode  degeneracy.  The  fit  for 
f±  =  1.2374  x  10 12  dyne/ cm :  and  <r  =  0.2703  has  a  (%)  =  0.0124.  This  is  sufficient  to  determine  ft  to  about  0.01%  and  cr  to 
about  0.05%.  There  are  no  corrections  so  these  values  are  absolute. 


n 

/,  (MHz) 

L  (MHz) 

%  error 

0U) 

1 

0.775706 

0.775707 

-0.000138 

(6, 1).  (1. 1),  (4. 1). (4.2). (7. 1) 

6 

0.819567 

0.819983 

-0.050778 

(5.1).  (3.1).  (5. 2).  (8.1).  (2.1) 

11 

1.075664 

1.075399 

0.024614 

(1.2).  (7.  2).  (6. 2) 

14 

1.198616 

!.  198505 

0.009239 

(5. 3).  (2. 2).  (3. 2).  (8. 2).  (3. 3).  (8. 3).  (2.  3) 

21 

1.217375 

1.2;  7850 

-0.039042 

(1.3).  (6,  3).  (7.  3).  ( 1 . 4).  (6.  4).  (7. 4).  (4, 3) 

28 

1.440760 

1.440750 

'  0.000712 

(5.4) 

29 

1.527080 

1.526474 

0.039695 

(5. 5).  (8.  4).  (3. 4),  (5. 6).  (2. 4) 

34 

1.558358 

1.558848 

-0.031448 

(5. 7).  (5. 8).  (5.  y).  (3. 5).  (8. 5).  (2. 5).  (3. 6).  (8. 6).  (2. 6) 

43 

1.580067 

1.579871 

0.012426 

(6.  5).  (7.  5).  (7.  6),  ( 1.  5).  (4.  4).  ( 1 . 6),  (6,  6).  (4.  5 ).  (4. 6) 

Determination  of  high  temperature  moduli  to  support  the 
broad  DOE  initiative  in  heat  treatment  distortion 

5120  steel  -  ground  parallelepiped 

free  dimensions  are  dl,  d2,  d3  (initial=  0.41205,  0.33973,  0.26380)  cm 
using  10  order  polynomials  mass=  .2875  gm;  density=  7.785  gm/cc 


n 

f-expt 

f-calc 

%err 

C11  C44 

df/d(moduli) 

1 

.340010 

.339607 

-.12 

.00 

1.00 

2 

.447200 

.447109 

-.02 

.13 

.87 

3 

.490240 

.490366 

.03 

.16 

.84 

4 

.541860 

.541482 

-.07 

.02 

.98 

5 

.557460 

.556663 

-.14 

.00 

1.00 

6 

.591740 

.592306 

.10 

.05 

.95 

7 

.608670 

.608417 

-.04 

.31 

.69 

8 

.615590 

.615290 

-.05 

.02 

.98 

9 

.623100 

.622662 

-.07 

.07 

.93 

10 

.642700 

.642897 

.03 

.12 

.88 

11 

.675510 

.676384 

.13 

.10 

.90 

12 

.683560 

.683316 

-.04 

.05 

.95 

13 

.690650 

.690436 

-.03 

.16 

.84 

14 

.746720 

.747335 

.08 

.10 

.90 

15 

.753860 

.754125 

.04 

.05 

.95 

16 

.827430 

.827569 

.02 

.06 

.94 

17 

.848300 

.849153 

.10 

.03 

.97 

18 

.855320 

.854953 

-.04 

.08 

.92 

19 

.870540 

.870865 

.04 

.12 

.88 

20 

.878390 

.878891 

.06 

.11 

.89 

21 

.882080 

.881734 

-.04 

.33 

.67 

22 

.884870 

.884342 

-.06 

.21 

.79 

23 

.887380 

.887639 

.03 

.20 

.80 

24 

.891190 

.891892 

.08 

.04 

.96 

25 

.916680 

.916689 

.00 

.31 

.69 

rms  error*  .0680  %, 

Fitted  values  for  elastic  constants  and  dimensions: 

Young’s  Modulus  *  30.39x10*  psi  (209.55  GPa) 

Shear  Modulus  =  11.82x10*  psi  (81 .50  GPa)  Poisson’s  ratio*  0.285 
dl  =  0.4121  cm  d2  =  0.3395  cm  d3  =  0.2639  cm 

Perturbation  analysis  of  error  sensitivity  :  (goodness  of  fit) 

X2  increased  2%  by  the  following  %  changes  in  independent  parameters 
C11  C44  dl  d2  d3 

.48  .41  -.42  -.42  -.4 


=>  errors  on  fitted  values  ~  0.5% 
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nn,  np,  it, 

rho»  6, 

10,  1, 

6.946 

i 

f  ex 

fr 

%err 

wt 

k 

i 

11 

33  23 

12 

44 

66 

1 

0.5501  01 

0.548752 

-0.25 

0.0 

4 

1 

0.01 

0.01  0.00 

0.00 

0.17 

0.33 

2 

0.655544 

0.653127 

-0.37 

0.0 

4 

2 

0.01 

0.00  0.00 

0.00 

0.47 

0.02 

3 

0.767276 

0.767542 

0.03 

1.0 

3 

1 

0.02 

0.01  -0.01 

0.00 

0.00 

0.48 

4 

0.840904 

0.840137 

-0.09 

1.0 

1 

1 

0.35 

0.17  -0.18 

0.01 

0.01 

0.15 

5 

0.842859 

0.840724 

-0.25 

1.0 

8 

1 

0.02 

0.01  0.00 

0.00 

0.48 

0.00 

6 

0.870588 

0.872288 

0.20 

1.0 

2 

1 

0.02 

0.01  -0.01 

0.00 

0.48 

0.00 

7 

0.682324 

0.882338 

0.00 

1.0 

6 

1 

0.41 

0.03  -0.05 

-0.06 

0.17 

0.00 

8 

0.914348 

0.916268 

0.21 

1.0 

7 

1 

0.31 

0.27  -0.24 

0.02 

0.04 

0.10 

9 

0.962734 

0.985844 

0.32 

1.0 

6 

2 

0.35 

0.06  -0.11 

0.01 

0.19 

0.00 

10 

1.017951 

1.018350 

0.04 

1.0 

5 

1 

0.40 

0.45  -0.41 

0.06 

0.00 

0.00 

11 

1*.  031479 

1.028883 

-0.25 

1.0 

1 

2 

0.07 

0.30  -0.09 

0.00 

0.14 

0.09 

12 

1.064057 

1.062296 

-0.17 

1.0 

2 

2 

0.36 

0.04  -0.09 

-0.01 

0.19 

0.01 

13 

1.070964 

1.073351 

0.04 

1.0 

7 

2 

0.10 

0.19  -0.04 

-0.01 

0.12 

0.14 

14 

1.071533 

1.072243 

0.07 

1.0 

5 

2 

0.47 

0.31  -0.27 

-0.01 

0.00 

0.00 

15 

1.091915 

1.091434 

-0.04 

1.0 

3 

2 

0.09 

0.45  -0.21 

0.02 

0.04 

0.12 

16 

1.101242 

1.101298 

0.01 

1.0 

2 

3 

0.05 

0.02  -0.02 

0,00 

0.21 

0.24 

17 

1.119834 

1.116659 

-0.28 

1.0 

8 

2 

0.37 

0.05  -0.08 

-0.02 

0.18 

0.01 

18 

0.000000 

1.135942 

0.00 

0.0 

5 

3 

0.63 

0.04  -0.02 

-0.14 

0.00 

0.00 

19 

1.363975 

1.164571 

0.05 

1.0 

4 

3 

0.02 

0.01  0.00 

0.00 

0.34 

0.13 

20 

1.186626 

1.188038 

0.10 

1.0 

8 

3 

0.08 

0.03  -0.04 

0.00 

0.19 

0.23 

23 

1.237561 

1.235048 

-0.20 

1.0 

5 

4 

0.32 

0.38  -0.26 

0.05 

0.00 

0.01 

22 

1.264551 

1.261810 

-0.22 

1.0 

3 

3 

0.07 

0.06  -0.04 

-0.01 

0.40 

0.01 

23 

1.296307 

1.294200 

-0.16 

1.0 

1 

3 

0.26 

0.33  -0.20 

0.01 

0.03 

0.08 

24 

1.317507 

3.316531 

-0.07 

1.0 

7 

3 

0.30 

0.31  -0.20 

0.00 

0.03 

0.06 

25 

1.321605 

1.320288 

-0.10 

1.0 

5 

5 

0.46 

0.11  -0.02 

-0.07 

0.02 

0.00 

26 

1.329940 

1.326090 

-0.29 

1.0 

6 

3 

0.45 

0.11  -0.15 

-0.02 

0.11 

0.00 

27 

1.354530 

1.354440 

-0.01 

1.0 

1 

4 

0.26 

0.07  -0.08 

-0.02 

0.11 

0.16 

28 

1.357242 

1.358113 

0.06 

1.0 

7 

4 

0.24 

0.16  -0.15 

0.01 

0.15 

•0.09 

29 

1.410719 

1.410488 

-0.02 

1.0 

6 

4 

0.27 

0.09  -0.09 

-0.02 

0.24 

0.00 

30 

1.471757 

1.474156 

0.16 

1.0 

6 

5 

0.27 

0.15  -0.15 

0.01 

0.21 

0.01 

31 

1.510419 

1.508057 

-0.16 

1.0 

3 

4 

0.16 

0.02  -0.03 

-0.01 

0.03 

0.34 

32 

1.511018 

1.515919 

0.32 

1.0 

7 

5 

0.35 

0.14  -0.04 

-0.06 

0.04 

0.07 

33 

1.512835 

1.516014 

0.21 

1.0 

6 

6 

0.13 

0.04  -0.03 

-0.01 

0.27 

0.11 

34 

1.516509 

1.518927 

0.16 

1.0 

1 

5 

0.23 

0.20  -0.09 

-0.02 

0.14 

0.05 

35 

1.564802 

1.560320 

-0.29 

1.0 

1 

6 

0.11 

0.06  -0.03 

-0.01 

0.27 

0.10 

36 

1.566096 

1.567781 

0.11 

1.0 

5 

& 

0.23 

0.21  0.01 

0.02 

0.03 

0.00 

37 

3.586001 

1.585657 

-0.02 

1.0 

8 

4 

0.15 

0.09  -0.06 

-0.01 

0.31 

0.01 

38 

1.592957 

1.590255 

-0.17 

1.0 

4 

4 

0.14 

0.06  -0.06 
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orthorhombic  elastic  constants  (cl  1 ,c22 ,c33,c23,cl 3, cl 2,c44,c55,c66) : 

2.65932  2.65932  2.58209  0.99080  0.99080  0.63820 

0.67740  0.67740  0.58697 

dimensions 

0.11436  0.12371  0.10992 

mass  error  •  0.00000e+00  % 

rms*  .  0.16522  percent 

chisquere  increased  2%  by  %  change  in  x's  of 
■  -0.25  -0.14  -0.66  -1.70  0.01  -0.01  0.00  0.00  0.00 

chisquere  increased  2%  by  %  change  in  x's  of 
0.06  -0.31  -0.28  0.53  0.00  0.02  0.00  0.00  0.00 

chisquare  increased  2*  by  %  change  in  x's  of 
-0.12  0.03  -0.13  0.39  0.05  0.00  0.00  0.00  0.00 

xohmrq&l  ctss  time  62.350  seconds 

cpu*  66.355  i/o»  7.158  mem«  6.837 


all  done 


Elastic  constants  of  tantalum  at  room  temperature. 


T(  K) 

P  (g/cm3) 

C,1  (GPa) 

cn  (GPa) 

C44  (GPa) 

3001 

16.678 

266.7 

160.8 

82.5 

3001 

16.678 

266.8 

161.4 

82.5 

3002 

16.633 

260.9 

157.4 

81.8 

29 83 

16.626 

260.2 

154.5 

82.6 

2954 

16.641 

266.3 

160.5 

82.8 

^.I.  Bolef,  J.  Appl.  Phys.  33,  2311  (1962). 

2F.H.  Featherstone  and  J.R.  Neighbours,  Phys.  Rev.  130,  1324  (1963). 

3N.  Soga,  J.  Appl.  Phys.  37,  3416  (1966). 

4Euler  angles  (defined  in  Roe’s  convention)  were  determined  to  be  a  =  138.3°, 
0  =  29.7°,  and  j  =  155.1°  by  RUS  and  a  =  135°,  0  =  33°,  and  7  =  158°  by 
X-ray. 


Vibrations  of  a  Rectangular  Parallelipiped 


*Not  a  plane  wave  among  them 

‘Sufficiently  complex  to  provide  all  the  elastic  moduli 

‘Completely  understood  mechanics  problem 
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John  William  Strutt,  the  Baron  Rayleigh  tried  to  do  this 
computation.  Without  a  90MHz  Pentium,  he  found  that 


In  the 

case  of  a  short  rod  and  of  a  particle  situated  near  the  cylindrical 
boundary,  this  lateral  motion  would  be  comparable  in  magnitude 
.with  the  longitudinal  motion,  and  could  not  be  overlooked  without 
risk  of  considerable  error. 

226.  The  problem  of  a  rectangular  plate,  whose  edges  are 
free,  is  one  of  great  difficulty,  and  has  for  the  most  part  resisted 
attack1.  ”  *  4 


Even  with  a  Pentium,  if  you  try  this  using  finite  element  methods, 
the  computation  time  goes  like  the  cube  of  the  numerical  accuracy 
and  you  can’t  compute  as  well  as  you  can  measure  in  a 
reasonable  time  on  a  reasonable  computer.  However,  if  you  are 
careful,  and  smart,  as  were  Orson  Anderson  and  his  postdoc 
Harold  Demarest  at  Bell  labs  30  years  ago,  then . 


Computation  of  resonances 
from  Migliori  et.al.  Physica  B  183,1,1993 


The  procedure  for  solving  the  direct  problem 
for  an  arbitrarily  shaped  elastic  solid  with  vol¬ 
ume  V ,  elastic  tensor  cijkI,  density  p,  and  with  a 
free  surface  S  begins  with  the  Lagrangian 

L  =  f  (KE  -  PE)  dV  (4) 

V 

where  the  kinetic  energy,  KE,  is  given  by 
KE  =  |  pa> 2«, 2 ,  (5) 


and  the  potential  energy,  PE,  by 
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Following  Hamilton,  we  allow  m,  to  vary  arbit¬ 
rarily  in  the  volume  V  and  on  the  surface  S 
+  8m,)  and  calculate  the  variation  8 L  in 

L.  The  result  is 

8L  =  J  (left  side  of  eq.  (8)), 8m,-  dV 
v 

+  J  (left  side  of  eq.  (9)),&u,  dS  (7) 

s 

The  immediate  results  are  two  equations,  the 
elastic  wave  equation  and  the  vanishing  of 

surface  traction 

pct>2wf.  +  Cijkluk  lj  =  0  ,  (8) 


njciikiukj  =  0  W 

where  {w,}  is  the  unit  outer  normal  to  S. 

Because  of  the  arbitrariness  of  5m,  in  V  and  on 
S ,  the  m/s  which  correspond  to  stationary  points 
of  L  (i.e.  8 L  =  0)  must  satisfy  eq.  (8)  in  V  and 
eq.  (9)  on  5.  There  are  no  such  m/s,  of  course, 
unless  a)2  is  one  of  a  discrete  set  of  eigenvalues, 
the  normal  mode  frequencies  of  free  vibration  of 
the  system. 


Following  the  Rayleigh-Ritz  prescription,  we 
expand  the  displacement  vector  in  a  complete  set 
of  functions  {<£A}, 

U,  =  «A,^A  >  (10) 

and  choose  as  our  basis  functions  powers  of 
cartesian  coordinates: 

*a  =  x'ymz\  (11) 

where  A  =  (/,  ra,  n)  is  the  function  label,  a  set  of 
three  nonnegative  integers.  After  substituting 
eq.  (10)  into  eq.  (4),  we  obtain  (a  becomes  a 
column  vector) 


L  =  \<x)2aJEa  -  \aJra 


(12) 
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Potato  Sphere  Ellipsoid  Spheroid  Sphere  Egg  Herrii 


FIG.  8.  Frequency  spectra  of  a  number  of  objects  in  the  potato  family.  The 
seven  stations  correspond  to  shapes  as  labeled,  with  semiaxes  as  given  in 
Table  II.  The  sphere  frequencies  agree  well  with  those  in  the  literature  for 
these  material  parameters  (Poisson’s  ratio  =  1/4).20  The  dimensional  pa¬ 
rameters  d\  +,dx  are  interpolated  linearly  between  the  seven  sta¬ 

tions  here.  Several  interesting  features  invite  comment.  First,  the  potato  has 
no  degenerate  lines,  because  of  its  low  symmetry,  and  the  sphere,  converse¬ 
ly,  has  few  lines  that  are  nondegenerate.  The  ellipsoid  has  no  degeneracies, 
and  the  spheroid,  the  egg,  and  the  hemisphere  (all  being  rotationally  sym¬ 
metric)  do,  but  never  more  than  doubly  degenerate  lines.  Small  deviations 
from  the  sphere  in  the  egg  direction  do  not  change  any  of  the  frequencies  to 
first  order,  because  d}  +  increases  as  much  a s  </3  _  decreases,  compensating 
one  another  as  far  as  affecting  resonant  frequencies  is  concerned.  As  in  sev¬ 
eral  other  figures,  apparent  avoided  crossings  on  this  plot  should  be  viewed 
with  suspicion  because  the  plotting  program  does  not  interchange  line  iden¬ 
tities  when  physically  the  modes  do,  in  fact,  cross.  Spectra  are  computed  for 
241  abscissa  values  here  and  elsewhere,  which  sets  the  scale  on  which  avoid¬ 
ed  crossings  may  be  spurious. 


Fig.  3 


Outline  sketch  of  the  stepplng-aotor  driven  system  used  to 
rock  the  sample  to  enable  detection  of  all  the  modes. 
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Transducer  construction 


Lithium  Niobate  Disc 


L 


Silvered  Kapton  Film 


Diamond 

Cylinder 


\ 


1.5  mm 


Silvered 
Kapton 
Electrical 
Lead 


Fig.  3.  Shown  is  a  schematic  of  the  diamond/polyimide/ 
LiNbO,  composite  transducer  used  for  all  the  measurements. 


RUS  can  be  used  at  temperatures  as  high 
as  1800  C  (O.  Anderson  et.  al.)-  For  more 
moderate  temperatures  (700  C),  both 
frequency  and  attenuation  data  can  be 
acquired  using  all-metal  diffusion  bonded 
LiNb03/diamond  transducers 


HighTemperature  RUS  System 

Piezoelectric  material :  LiNb03 :  Tc~1 197,  TrrM260 


Transducers 


0.001"  thick  SS  mesh 


Backload7 


l 


Diffusion 
Metal  bonds 


Probe  assembly 
for  furnace 


Parallelepiped  sample  held 
lightly  between  transducers 


Stainless  to  stainless  spot  weld 


Alumina  ceramic 


Mesh 


Stainless  steel 


Figure  7 


Silver  diffusion  bonded  0.375"O  PZT-5A 
Transducer  with  Alumina  wear  plate  and 

alumina  backload 


The  bare  transducer  has  a  2MHz 
fundamental  compressional  mode 
This  technology  won  a  1993  RD100  award 
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RUS  Cell  with  Silver  diffusion  bonded 
transducers,  using  LiNb03  and  an  alumina 
cylinder  backload.  This  cell  can  reach  500C 


The  transducers  are  mounted  beneath  the 
0.500"0  304SS  screens  that  the  sample 
touches.  Only  the  bottom  one  is  visible. 


Amplitude  (art),  units) 


Resonances  of  a  51 20  steel  RPR  at  38C  and  378C 
measure  using  metal  diffusion  bonded  LiNb03/Alumina 
transducers 
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The  mechanical  resonances  of  a  solid  are  determined  by 
the  bulk  thermodynamic  properties  as  well  as  the 

morphology. 


NdjFe,^  Magnet  material  - 

resonance  frequencies  near  the  Curie  temperature 


Temperature  °C 

Hard  magnets  are  hard  because  impurities  and  grain 
boundaries  in  very  high  permeability  ferromagnets  affect  the 
ability  of  magnetic  domains  to  track  an  applied  field.  If  the 
domains  can’t  move,  the  magnet  is  hard.  The  entire  trick  with 
Nd2Fe14B  was  to  determine  the  changes  to  the  intergranular 
region  that  stopped  the  domains  from  moving. 

RUS  can  extract  useful  information  relating  to  the  effect  of 
the  intergranular  composition  on  the  overall  magnetic 
properties  of  Nd2Fe14B. 
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Use  of  RUS  to  determine  the  spin  reorientation 
temperature  of  Nd2Fe14B  obtained  from  GM  Research 
Laboratories 


Temperature  (K) 

FIG.  1.  (a)  Longitudinal  and  (b)  transverse  elastic  constants  of  Pr2Fc,4  B  and 
Nd2FeuB  in  the  100-200  K  temperature  range. 


The  phenomenon  of  spin  reorientation,  depicted  above, 
can  occur  at  technologically  inconvenient  temperatures, 
changing  the  important  properties  of  a  permanent  magnet 
used  under  extreme  condition  (such  as  under  the  hood  of 
an  automobile).  RUS  provides  a  convenient  probe  for  such 
effects  that  is  sensitive  to  the  average  concentration. 
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Resonant  ultrasound  studies 
-superconductivity 


One  of  the  clearest  examples  ever 
of  the  discontinuity  in  elastic  properties 
at  a  superconducting  phase  transition 


663.5 
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663.3 


YNi2B2C  Single  Crystal  Flake  AL227 
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Resonant  ultrasound  studies 
-valence  transition 

Valence  transition  at  40K  and 
phase  transition  of  unknown 
origin  at  66K 
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U2Zn17 

Rhombohedral,  macroscopic 
hexagonal  symmetry 

Mirror  twins,  stacking  faults 

Kondo  system  with  enhanced 
specific  heat  and  T3term 


Fig.  II.  Spcuhc  heal  minus  the  lattice  contribution  o/LjZn,.  plotted  afJlT  versus  Ta.  Data  are  from 
Ott  ci  I  i  VK4bL  The  broken  lines  indicate  the  temperature  dependences  above  and  below  (ha  phase 


IIKW  VLLECTKON  ACTINIDE  MATERIALS 


Commensurate  Antiferromagnet 

Nuclear  and  magnetic  unit  cells 

the  same  3 

« 

Easy  axis  in  hexagonal  plane 


mognelfC  susceputahty 

KX)  u,Znlf  / 


Tn=9.8K 


£  60  F  /  2sL©*X<* 

■  K 

X  •  / 

40 

/>  20 


150  200  250  300 
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Fig.  25.  Temperature  dependence  t»f  ihe  invert  magnetic  susceptibility  /  1  of  UjZn,,  between  1.3  and 
HA)  k  The  inset  sho* »  /}  / 1  on  an  cvicmicd  temperature  scale  below  30  K..  (Data  are  taken  from  Oil  et 


Fig.  51.  Temperature  dependence  of  the  electrical  resistivity  of  UjZn,-  between  1.2  and  300  K.  Note  the 
different  scales  for  different  temperature  intervals.  (Data  are  from  On  et  al.  I984b.| 
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Antiferromagnetic  Ordering  of  U 


U2Zn,7 


Fig.  113.  Proposed  magnetically-ordered 
structure  of  V2Znn  on  the  basis  of  neutron- 
diffraction  results  by  Cox  et  aL  (1986^  Open 
and  closed  circles  denote  atoms  with  oppositely 
oriented  moments.  The  moment  directions  are 
wtthin  (he  {uual  plane. 


Stacking  Faults 
Twins 

Remnants  of  high  T  hexagonal  structure  in  low  T  rhombohedral 
phase 

The  above  can  lead  to  weak  ferromagnetism  and  structural  effects 
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The  mechanical  resonances  of  a  solid  are  determined  by 
the  bulk  thermodynamic  properties  as  well  as  the 

morphology. 
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Hard  magnets  are  hard  because  impurities  and  grain 
boundaries  in  very  high  permeability  ferromagnets  affect  the 
ability  of  magnetic  domains  to  track  an  applied  field.  If  the 
domains  can’t  move,  the  magnet  is  hard.  The  entire  trick  with 
Nd2Fe14B  was  to  determine  the  changes  to  the  intergranular 
region  that  stopped  the  domains  from  moving. 

RUS  can  extract  useful  information  relating  to  the  effect  of 
the  intergranular  composition  on  the  overall  magnetic 
properties  of  Nd2Fe14B. 


NdjFe^B  Magnet  material  - 

resonance  frequencies  near  the  Curie  temperature 
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Basic  Research  in  High-Temperature 
Superconductors  Leads  to  a  Unique  and 


Particle  Accelerators 


*An  insoluble  problem: 

*Reason:  it  is  just  like  chess  but  without  the  little  men 

*Proof:  No  two  particle  accelerators  are  alike 

No  minimalist  approach  is  possible  because  the  problem 
is  too  il-defined  and  complex.  But  nevertheless,  great 
physics  comes  of  it,  at,  of  course,  great  expense.  One 
cannot  give  up  the  physics,  one  can  only  lament  the 
expense 


The  Josephson  Voltmeter 

*Benchtop  physics  with  a  few  basic  variants 

*Fantastic  macroscopic  manifestation  of  quantum 
mechanics,  as  wonderful  and  important  as  any  bit  of 
physics 

*Cheap  (costs  less  than  a  car,  is  smaller  than  a  Kitchen) 
*Redefines  the  standard  for  Voltage  (i.e.  it  is  very  useful) 


/* 

/ 


/ 


112 


The  Crescent  Wrench: 

*A  problem  that  can  be  solved: 

*Reason:  The  problem  is  simple,  as  are  the  approaches 
for  solution 

*Proof:  There  is  only  one  basic  design  left  at  the  hardware 
store 


Chess: 

*An  insoluble  problem: 

*Reason:  it  is  too  complicated  and  messy,  with  too  many 
possibilities 


*Proof:  People  still  play  chess 
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Elastic  moduli  of  U2Zn17 

Rectangular  parallelipiped,  crystal  axes  not  parallel  to  faces  but 
determined  by  RUS 


c33 

ell 

c23 

c44 

c66 


Data  taken  in  two  separate  runs.  Failure  to  get  perfect  match  is 
typical  if  morphological  flaws  exist. 

Shear  moduli  ±.5%,  diagonal  moduli  ±2%,  off  diagonal  moduli  ±4% 


c44  is  unaffected  by  Neel  transition,  all  other  moduli  soften 
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Resonant  frequency  of  single-crystal  fragment 
Undetermined  mix  of  cM  and  q* 

Undetermined  orientation  of  field  with  crystal  axes 


Variation  of  Neel  temperature  and  softening  with  field  and  with 
c^/Cee  ratio  (undetermined)  on  several  samples  with  different  mixes 

Tn  essentially  independent  of  field  direction,  suggesting  free 
energy  couples  as  M2H2 


Effects  above  the  Neel  transition  in  applied  field 
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The  above  is  indicative  of  a  weak  phase  transition  of  unknown 
origin 


This  effect  disappears  above  about  70K 


Frequency(kHz) 
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UCu5  05  annealed  polycrystal  showing 

the  clear  slope  break  in  the  shear  modulus 
at  the  antiferromagnetic  transition 
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Simple  model  for  the  unusual  band-edge 
density  of  electronic  states  in  FeSi 


Temperature  (K) 


FIG.  1.  Magnetic  susceptibility  of  FeSi.  Open  circles:  experi¬ 
mental  points  after  Jaccarino  et  al.  (Ref.  5).  A  low-temperature 
Curie  tail  was  subtracted  from  the  data  as  described  in  Ref.  5. 
Solid  line:  calculation  using  the  model  density  of  states  shown 
in  the  inset  with  parameters  £g=950  K,  W =750  K,  and 
g=4.40  states/cell.  Dashed  line:  calculation  using  parameters 
Eg  =  1 100  K,  W = 500  K,  and  g = 4. 20  states/cell. 
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Electronic  Structure  of  the  Narrow-Gap 
Semiconductor  FeSi  using  RUS 


In  order  to  fit  the  data  better,  we  consider  a  deforma¬ 
tion  potential  coupling  which  explicitly  includes  the  con¬ 
tribution  of  conduction  electrons  to  the  elastic  moduli 
through  a  rigid  two-band  model  (E(k)  =  E°(k)  -I-  dr(k)er, 
where  dr(k )  is  defined  as  dE(k)/dsr>  and  is  a  symmetry 
strain) 


Consider  the  free  energy  for  conduc¬ 
tion  electrons  with  band  index  i  and  energy  El(k): 


kBT^\n 

i .  k 


1  4-  exp 


where  (x  is  the  chemical  potential.  Explicitly  calculating 
the  symmetry  elastic  moduli,  cr  =  d2F/dE2 ,  and  assum¬ 
ing  conservation  of  the  total  number  of  quasiparticles  [5] 
yields 

-A- £</?<*)/*< l  -ft) 

Kb1  k 

1  (I  kdr(k)ft(\-ft))2 

+  k»T  ItMl-ft) 


(3) 


Thermal  expansion  of  the  semi-metal  CoSi 
provides  background  thermodynamic 
information  to  enable  extraction  of  the 
electronic  effects  in  FeSi 


Temperature  (K) 

Fig.  1.  Linear  expansion  coefficient,  a,  versus  temperature  for 
Fei-jCo^Si  with  x  =  0,  0.02,  0.05,  and  1. 
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Elastic  Constants  of  FeSi 
Physica  B  478,199  (1994) 


Fig.  1.  Elastic  moduli  of  FeSi  as  a  function  of  temperature.  The 
dashed  curves  are  fits  to  the  Varshni  function  (Eq.  (1);  for 
cu,s=  70.7  GPa,  r  =  365  K;  c44,  s  =  39.5  GPa,  t  =  366  K; 
and  for  l/2(cu  -  c12),  s  =  40.0  GPa,  r  =  374  K).  The  solid 
curves  are  fits  using  a  deformation  potential  coupling  model. 
The  parameters  for  these  fits  are  given  in  Table  1. 


1 000/Q 


Ultrasonic  Attenuation  of  FeSi 
Physica  B  478,199  (1994) 


Modulus 

d(K) 
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17 

4951 
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7 

10600 
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Fig.  2.  1  /Q,  which  is  proportional  to  attenuation,  is  plotted  for 
one  of  the  measured  resonance  frequencies.  The  peak  in  attenu¬ 
ation  occurs  at  the  activation  energy  (d/2). 


Structural  phase  transition  in  La2.xSrxCu04 

RUS  simultaneously  sees  collapse  of  c66 
and  no  effect  on  cA4  in  a  2mm  single  crystal 


Temperature(K) 


Temperature(K) 


La2Cu  Q* 
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Numerical  analysis  of  the  motion  (6)  establishes  that  the  eigenmode  shown  in  Fig.  1 
depends  almost  purely  on  c^.  Absent  dynamical  effects,  we  can  treat  the  temperature 

dependence  of  c^  with  a  conventional  Ginsburg-Landau  (G-L)  Hamiltonian.  A 

simplified  version  of  such  a  Hamiltonian,  containing  all  the  essential  physics  is 

F  =  1/2  a(T-Tc)  q2  +  1/4  P  q^  +  1/2  c ^  e^2  +  8  e^  q2  ( 1 ) 

where  e^  is  the  shear  strain,  q  is  one  component  of  the  order  parameter,  and  a,  P,  and 

8  are  parameters.  The  last  term  in  Eq.  1  is  the  quadratic  coupling  between  strain  and 
order  parameter.  Minimization  of  this  expression  with  respect  to  q  determines  q(T). 
Using  q(T)  we  can  then  calculate  c^^(T)  by  taking  the  second  derivative  of  F  with 

respect  to  strain  to  obtain 

c66(T)  =  c66  forT>Tc 

(2) 

c66^  =  c66-252/P  forT<^c’ 

where  q  is  zero  above  T  and  non-zero  below.  Eqs.  2  are  the  usual  G-L  result, 

v 

discussed  more  completely  by  Rehwald  (7).  We  see,  therefore,  that  simple  quadratic 
coupling  and  no  dynamics  produces  only  a  step  discontinuity  in  c^  at  the  structural 

phase  transition.  This  is  not  what  the  data  show.  The  data  fit  a  Curie-Weiss  (C-W) 
softening  of  the  form 

WT>  =  <W1-S/<T-Tc>>  (3) 

where  Tfi  is  223K  and  the  fit,  shown  as  the  solid  line  in  Fig.  1,  is  accurate  to  0.2%  over 

more  than  a  decade  in  g/(T-T  ). 

c 

Gaussian  fluctuations  of  the  order  parameter  (8),  self-consistent  phonons  (9),  and 
linear  coupling  between  strain  and  order  parameter  (7)  all  yield  a  C-W  dependence  for 
c66'  ^or  Gaussian  fluctuations,  the  critical  exponent  for  the  specific  heat  and  for  the 

elastic  moduli  is  X  =  2  -  d/2,  where  d  is  the  dimension  of  the  order  parameter.  In  our 
system,  the  order  parameter  is  two  dimensional,  thus  the  critical  exponent  (the 
exponent  of  1/(T-T  ) )  is  unity,  in  agreement  with  the  data.  However,  our  C-W  fit  is 

V 

over  a  temperature  range  of  about  80K  (g=1.47K).  This  is  a  very  large  range  for 
fluctuations  to  be  important.  We  can  make  an  upper-bound  estimate  for  the  fluctuation 
regime  (8)  by  using  a  few  lattice  spacings  for  the  coherence  length,  and  by  using  a  1% 
(SrTiO^  has  about  a  10%  modulus  discontinuity  at  its  phase  transition  temperature) 

modulus  discontinuity  to  obtain  a  fluctuation  range  of  about  IK,  comparable  to  the 
region  in  Fig.  2  where  the  ultrasonic  attenuation  increases  sharply.  Thus  it  appears 
very  unlikely  that  2-D  Gaussian  fluctuations  can  explain  what  we  observe. 


zone-edee  soft  o  l  !  anhannonic  potential  associated  with  the 

For  rhiccn^  f VT  d  f  UlC  °  octahedra  can  Produce  C-W  modulus  softening  (9) 
For  this  son  of  treatment  to  be  successful,  the  zone-edge  soft  mode  must  be  linearlv 

bv^nsmin^6  ZOnKe'center  acoustic  mode-  The  1-D  treatment  in  ref.  (9)  deals  with  dm 
?ond!f  ?Su  anharmonic  sPnng>  used  in  the  shell-model  to  develop  the  self- 

conShmi/t^011  £hspefs!on  curve>  in  series  with  the  ion  cores.  Thus  this  spring 
contnbutes  to  the  potential  energy  for  any  value  of  k,  the  phonon  wave  vector.  § 

hlanoh  °  SCatlerinf-  measurements  (3)  show  that  the  soft  mode  is  pan  of  the  phonon 
branch  corresponding  to  c^,  not  c66.  Without  some  linear  coupling  term  to  the  c„ 

not.easy  t0  see  the  applicability  of  self-consistent  phonons  to 
^U  U s°f[ening  of  concern  here.  Were  such  a  term  to  be  added,  the 

enera  W°U  1  t0  b?come  explicitly  3-D,  and  because  both  the  coupling  and  the 

energies  would  depend  on  the  anharmonicity,  the  C-W  exponent  would  likely  be  lost 

The  third  possibility  we  consider  here  is  the  replacement  of  quadratic  coupling  with 
inear  coupling  (  for  T>Tc,  the  inclusion  of  the  quadratic  term  has  no  effect  with  or 

wnhout  the  linear  term  present)  between  order  parameter  and  strain  in  Eq.  1.  This 
F  =  1/2  a(T-Tc)  q2  +  l/4  p  q4  +  i/2  c66  e62  +  y  e6  q  (4) 


c66(T)  =  c66  -  ^(“(T-T- 


for  T>T 


2  O) 

C66(T>  =  c66  '  r/(2a(T-T))  forT<Tc, 

as  required  to  fit  the  data  of  Fig.  1.  To  justify  a  linear  coupling  term,  the 
a1.86*r.14Lu04  CT^stal  must either  non-linear  or  not  tetragonal. 

In  Fig.  4  we  plot  some  of  the  lowest  eigenfirequencies  of  the  Laj  g6Sr  14Cu04  crystal 

5  P1011^  ****  Iowest  two  eigenmodes  on  an  expanded  scale, 

2"  Cr0SS1T  0f  SCVCral  pcrcent-  Note  “  Fig-  4,  many  avoided 

th^r^  TT  The|.mechanical  Lagrangian  for  analysis  of  the  resonances  of 

h!~  I$  ba$ed  T  a  lmfar  tetrag°nal  sobd  d.6).  This  model  produces  eight 

2?^^  SSCS  °r  thC  modes-  ^  dimensions  of  the  crystal  are  fuch 

that  d0  CT0SS  as  temperature  is  varied  are  all  orthogonal  so  that  none 

crvshtalTnnndiCT0SSingS  Sh°Ul<1  OCCur*  Their  existence  can  explained  only  if  the 
crystal  is  non-lmear,  not  tetragonal,  or  has  excessive  preparation  errors. 
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Macroscopic  Ginsburg-Landau  picture 
Free  energy  expansion 

F  =  Fop  +  Fe  +  Fop-e 

F0p  =  a(Qj2  +  Q22)  +  u(Qj2  +  Q22)2  +  v(Qj4  +  Q24), 

where  and  Q2  are  the  order  parameters  corresponding 
to  the  two  directions  of  octahedral  tilt. 

Fe=  l/2c11(e12+e22)  +  c12eie2  +  c13(e1  +e2)e3 

+  1/2  C33e32  +  1/2  c^Ce^  +  e52)  +  1/2  c66e62 

where  cy  are  the  elastic  constants  and  ej  the  macroscopic 
strains. 

^op-e  =  ^xx^el+  e2^  +  ^zze3^Ql^  +  Q2^ 

+  kxye6^1^  “  Q2^)» 

where  are  the  symmetry-allowed  coefficients  from  the 
general  expansion 

F=  ^a,b,i  ^abeabQi^* 


Ginzburg-Landau  description  of  second  order  phase  transitions 
Linear  coupling  between  strain  e4  and  order  parameter  q 


F  =  J«(T  -  Tc)q2  +  ipc,4  +  lCi4e2  +  ye4g 

1 .  Minimize  F  with  respect  to  q. 

dF 

—  =  a(T-Tc)q  +  fiq3  +ye4  =0- 
dq 

Note  that  at  zero  strain,  7>7C,  qr=0,  while  for  7<7C  q#0. 

This  is  the  constraint  that  we  use  when  computing  the  elastic  response. 

2.  Compute  the  elastic  response  under  the  constraint. 
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c“<T>-0l-- 

Here’s  how!  Use  the  constraint  that  F  be  a  minimum  and  compute: 


^-[<x(7-7c)<7  +  |3<r>+Ye4]  =  0 


de4  a(T-Tc)  +  3fiqd 


d2q  _  3(3yg _ dq  _  -3pqr|  dq  j 

5e4  (a(T-Tc)  +  3$q2f  de4  y 
Now  compute  cu(T) 

g  =  (a(7-7c)  +  3Pg2  (jg  +  (a(T  -  7C),  +  P93)  g  +  c„  ♦  2y  g  +  re4 

Noting  that  if  7>TC  then  q=Q  and  that  we  are  evaluating  at  e4=0  we  get: 

dq  =  -y 
de4  a(T-Tc) 


giving 

c44(r)  =  C44  - 


«(r-rc) 


if  T>T, 


while  for  T<TC  we  find  that 


and  g2  _  ~a0~ ~TC) ,  producing  a  so-called 


0,4(7)  =  0,4  +  -  7I—  and  g*  =  ,  prc 

44  v  '  44  2oc(7-7c)  w  P 

Curie-Weiss  behavior  on  either  side  of  the  transition. 
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Let’s  now  look  at  quadratic  coupling  between  strain  and  order  parameter 

F  =  |a(7-  Tc)q2  +  |p<74  +  lcuej  +  ye4q2 

Minimizing  F  with  respect  to  the  order  parameter  q  yields 

=  (a(T  ~Tc)  +  2ye4)q  +  P<73  =  0 

which  has  as  solutions 

<7=0  above  Tc,  independent  of  e4 

and 

a2  -  +  for  T<TC  as  ling  as  the  strain  is  small. 

P 

The  effect  of  this  is  that  for 
T  >TC  c44  (7)  =  C44 

and  for 

Y2 

T<TC  c44(T)=cm-1~- 


Thus  a  step  discontinuity  occurs! 
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What  happens  when  coupling  is  turned  on; 


F  =  ^a(T-Tc)q2  +  lfiq,  +±C„eJ  +ye4<7” 

^-  =  0  =  04464  +  79" 

5e6 


T^TC  + 


ac^ 


if  n  =  1 


2y2 

=  if  n  =  2 

C44 

Thus  the  effect  is  to  either  renormalize  Tc  if  n=1 ,  or  if  n=2,  the  transition  can  become 
first  order.  Thus  elastic  coupling  can  have  either  weak  or  strong  effects  on  the 
underlying  physics.  Note  that  a  really  careful  treatment  will  always  produce  a  first  order 
transition,  albeit  weak,  if  n=2. 
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Now  to  make  predictions, 

minimize  free  energy  with  respect  to  strain  (the  e^'s  )  and  the 
order  parameters  (the  Qj's). 

This  minimization  predicts: 

1.  A  C44  =0 

2.  A  C}  1=  A  Cj2»  so  A  (cj 

3-  A  c66  »  A  C33  >Ac13>Ac11. 
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A .  Migliori  et  al.  /  On  techniques  for  measurement  of  the  elastic  moduli  of  solids 


*  Cll 

*  Cl  2 

*  C44 


Temperature  (K) 

Fig.  10.  The  three  elastic  moduli  of  a  single  crystal  of  SrTiO,  near  the  structural  phase  transition  are  shown  as  a  function  of 
temperature.  These  data  were  obtained  using  RUS. 
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Structural  phase  transition  in  a  1 .5mm 
single  crystal  of  pure  La2Cu04  at  527K 
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terials.  Deuterium  was  loaded  fron  the  gas  phase  in  the 
usual  way.8,16  Rectangular  parallelepipeds,  ~2  mm  on 
an  edge,  were  cut  from  larger  samples  using  an  electric 
discharge  machine.  Vacuum  fusion  analysis  of  one  of  the 
YDX  crystals  showed  x  —0. 10,  and  260  ppm  of  N  and 
2890  ppm  of  O  impurities.  Whereas  qualitative  results 
were  obtained  on  two  different  samples,  the  data  reported 
below  were  obtained  on  a  single,  rectangular- 
parallelepiped  sample  of  room-temperature  dimensions 
1.82X1.92  mm  X  2. 06  mm.  The  long  axis  of  the  paral¬ 
lelepiped  was  parallel  to  the  c  axis  to  within  0.5*,  as 
determined  by  x-ray  diffraction.  Since  a  crystal  of  hexag¬ 
onal  symmetry  is  elastically  isotropic  in  the  basal  plane, 
the  orientations  of  the  other  crystalline  axis  were  not 
determined  relative  to  the  parallelepiped  axes. 

Resonant  ultrasound  spectroscopy17  was  used  to  mea¬ 
sure  the  ultrasonic-attenuation  and  elastic  constants. 
With  this  technique  a  single-crystal,  rectangular- 
parallelepiped  sample  is  placed  comer-to-comer  between 
two  piezoelectric  transducers,  one  transducer  being  used 
for  generation  and  the  other  for  detection  of  ultrasonic 
vibrations.  By  sweeping  the  excitation  frequency,  a  large 
number  of  the  lowest  vibrational  eigenfrequencies  can  be 
investigated.  The  eigenfrequencies  are  related  to  the  elas¬ 
tic  constants,  while  the  inverse  of  the  Q  of  the  resonances 
is  related  to  the  vibrational  energy  loss  in  the  sample. 
Using  the  method  of  Visscher  et  al., 14  we  have  identified 
the  lowest  31  eigenfrequencies  of  one  sample  by  fitting 
the  measured  frequencies  to  a  set  of  elastic  constants. 
This  identification  is  needed  for  a  microscopic  interpreta¬ 
tion  of  the  results  as  discussed  below. 

III.  RESULTS  AND  DISCUSSION 

Figure  1(a)  gives  results  for  1  /Q,  for  a  single-crystal 
YD0  10  sample,  obtained  at  a  frequency  of  0.81  MHz, 
over  the  temperature  range  of  15-325  K.  Three  general 
features  are  apparent.  The  loss  increases  at  the  highest 
temperatures;  this  increase  may  be  associated  with  the 
long-range  diffusion  process,  or  other  processes  discussed 
below.  There  is  a  small  attenuation  effect  near  160  K,  in 
the  temperature  range  where  resistance  anomalies  have 
been  observed  and  associated  with  hydrogen  ordering.19 
The  attenuation  anomalies  in  this  temperature  range 
were  found  to  depend  on  the  rate  and  direction  of  tem¬ 
perature  change  and  have  not  been  well  characterized. 
The  main  effect,  of  course,  is  the  large,  broad,  asym¬ 
metric  attenuation  peak  with  the  maximum  near  87  K. 
The  solid  line  is  a  theoretical  fit  to  the  data  to  be  de¬ 
scribed  below.  Figure  Kb)  gives  data  similar  to  that  of 
Fig.  1(a),  but  for  a  different  vibrational  eigenmode,  the 
difference  to  be  discussed  below.  The  major  difference 
between  the  data  sets  for  the  two  eigenmodes  is  in  the 
magnitude  of  the  low-temperature  peak.  Figure  2  gives 
additional  data  at  other  frequencies  for  the  low- 
temperature  peak.  A  temperature-independent  back¬ 
ground  (1  /Q)bg  was  subtracted  from  the  two  highest- 
frequency  curves  before  plotting  the  data  in  Fig.  2. 
(\/Q)h  =0.00007  and  0.0001  for  the  1.4-  and  2.7-MHz 
curves,  respectively.  Since  the  peak  for  the  lowest  fre¬ 
quency  was  much  higher  relative  to  the  background,  no 


FIG.  1.  (a)  Ultrasonic  loss,  1/Q  vs  temperature  for  YD0  10  at 
a  frequency  of  0.81  MHz.  The  solid  line  represents  a  theoretical 
fit  to  the  data  described  in  the  text.  This  vibrational  eigenmode 
depends  almost  purely  on  Cu.  (b)  Ultrasonic  loss.  1/Q,  vs  tem¬ 
perature  for  YD0  ,0  at  a  frequency  of  0.83  MHz.  This  vibration¬ 
al  eigenmode  depends  almost  purely  on  Cw. 


FIG.  2.  Ultrasonic  loss,  1/Q.  vs  temperature  at  three 
different  frequencies.  The  solid  lines  represent  theoretical  fits  to 
the  data  described  in  the  text. 
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modes  by  fitting  our  measured  frequencies  to  a  set  of  elas¬ 
tic  constants  using  the  method  described  by  Visscher 
et  a/.24  In  general,  each  of  the  vibrational  eigenmodes 
depends  on  a  combination  of  elastic  constants.  However, 
the  lowest  eigenfrequency  and  a  few  of  the  higher  eigen- 
frequencies  depend  on  only  one  or  two  elastic  constants. 
This  fact  enables  us  to  obtain  a  rather  accurate  descrip¬ 
tion  of  the  temperature  dependence  of  some  of  the  elastic 
constants  by  simply  measuring  the  temperature  depen¬ 
dence  of  a  selected  number  of  eigenfrequencies.  In  addi¬ 
tion,  the  ultrasonic  attenuation  associated  with  such 
eigenfrequencies  is  obtained  from  these  measurements. 
In  the  present  study  the  emphasis  is  on  the  temperature 
dependence  of  the  ultrasonic  loss,  which  is  obtained  from 
the  inverse  of  the  Q  of  the  eigenmodes.  Those  modes 
which  depend  on  only  one  or  two  elastic  constants  are 
particularly  useful  for  making  a  microscopic  interpreta¬ 
tion  of  the  phenomena  responsible  for  the  ultrasonic  loss. 


Second,  the  loss  increases  as  the  temperature  approaches 
room  temperature;  this  increase  is  probably  associated 
with  the  long-range  diffusion  process.  Finally,  the  dom¬ 
inant  effect  is  the  broad  loss  peak  with  the  maximum  near 
50  K.  These  results  will  be  discussed  in  more  detail 
below. 

Figure  4  shows  data  for  ScH0  25  at  three  different  fre¬ 
quencies.  The  solid  lines  are  theoretical  fits  to  the  data  to 
be  described  below.  Figure  5  shows  similar  data  for 
ScD0>18.  The  different  frequencies  correspond  to  different 
vibrational  eigenmodes.  It  is  seen  that  the  attenuation- 
peak  position  shifts  little,  if  any,  with  a  change  in  measur¬ 
ing  frequency.  The  attenuation  peak  for  ScD0  ,8  occurs 
at  a  considerably  higher  temperature  than  is  the  case  for 
ScH0. 25- 

Measurements  were  also  made  on  an  undoped  Sc  crys¬ 
tal.  There  is  no  low-temperature  peak  in  this  case,  only  a 
small,  almost  temperature-independent,  attenuation. 


III.  RESULTS 


IV.  DISCUSSION 


Figure  3  shows  results  for  the  ultrasonic  loss  1  /Q  for 
two  different  vibrational  eigenmodes  in  ScD0  18.  Our 
analysis  shows  that  the  two  modes  of  Fig.  3  depend  only 
on  the  shear  elastic  constants  C 44  and  C^;  that  is,  the 
derivative  of  the  computed  eigenfrequency  with  respect 
to  the  elastic  constants  is  zero  for  the  other  three  elastic 
constants.  Further,  the  lower-frequency  mode,  at  0.73 
MHz,  is  almost  pure  Cu ;  the  dependence  of  the  higher- 
frequency  mode,  at  1.01  MHz,  on  elastic  constants  is 
weighted  about  40%  C44  and  60%  C^.  There  are  three 
main  features  shown  in  Fig.  3.  First,  there  is  some  evi¬ 
dence  of  an  ultrasonic-attenuation  anomaly  near  170  K, 
in  the  temperature  range  where  resistance  anomalies  have 
been  observed  and  associated  with  hydrogen  ordering.25 
The  attenuation  anomalies  in  this  temperature  range 
were  found  to  depend  on  the  rate  and  direction  of  tem¬ 
perature  change  and  have  not  been  well  characterized. 
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FIG.  3.  Ultrasonic  loss  i/Q  for  two  different  vibrational 
eigenmodes  in  ScDq  jg.  The  mode  at  0.73  MHz  is  almost  pure 
C44,  whereas  the  mode  at  1.01  MHz  is  about  40%  Cu  and  60% 


We  describe  our  results  in  terms  of  the  theory  of  two- 
level  systems  (TLS’s).  An  interstitial  such  as  hydrogen, 
which  may  occupy  either  of  two  nearby  interstitial  sites, 
can  be  described  as  a  two-level  system  with  an  energy 
splitting  E  ={Ej-+  A2)w1,  where  ET  is  the  tunnel  split¬ 
ting  and  A  is  the  difference  in  energy  of  the  two  wells 
(asymmetry).  This  formalism  may  be  used  even  cases 
where  tunneling  is  not  a  factor,  in  which  case  we  simply 
have  £r  = 0.  The  ultrasonic  loss  due  to  relaxation  is 
given  by26 

1  _  n°Dl  -cch2  f  E  }_cor_ 

Q  4Ck„T  *  \2kBT\  1+uV  ’ 

where  nQ  is  the  concentration  of  TLS’s,  D  =  bE  /bt  is  the 
variation  of  the  energy-level  splitting  with  respect  to  the 
ultrasonic  strain  e,  C  is  an  elastic  constant,  cj/Itt  is  the 
ultrasonic  frequency,  and  r  is  the  relaxation  time.  It  is 
usually  the  case27  that  6/1  /8e  »8£r/6e,  so  that 
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FIG.  4.  Ultrasonic  loss  1  /Q  for  three  different  frequencies  in 
ScHo  25.  The  solid  lines  are  fits  to  the  data  described  in  the  text. 
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RUS  and  non-destructive  testing-or  how  to  make  everyone  happy 
with  good  physics! 

Quality  control  and  the  means  to  test  for  it  have  become  of 
increasingly  greater  concern  to  American  industry,  and  therefore,  to 
government  and  other  funding  sources. 


Not  like  this 


141 


Detection  of  microcracking  in  a 
single  crystal  sample  of  EuB6 

A)  Resonances  of  the  pristine  sample 

B)  Resonances  after  a  small  stress  has  induced 

microcracking 
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Resonant  Ultrasound  Spectroscopy 
-nondestructive  testing 

Sphericity  error  shown  by  splitting 
of  degenerate  spherical  resonance 
in  a  Si3N4  ball  bearing  from  NASA 


Commercial  RUS  system  for  production  testing  of  automotive  oxygen  sensors. 
This  system,  as  well  as  other  RUS  systems,  is  built  and  marketed  by 

Quatrosonics  Corp.  Of  Albuquerque,  NM 


Oxygen  sensor  testing-flaw  selection  sensitivity 
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Production  testing  of  zirconia  automotive  oxygen  sensors  using 
high-speed  automated  resonant  ultrasound  spectroscopy  (RUS) 


146 


O  U 

m  ® 
CD  ^ 


in  o 

CM  O 

•  m 

o  o 
»n 

v*  7Z 


E 

o 

it 


0)  < 

©> 
CO  < 


o  ®s 
©  ^  O 

c  tc  l — 
CD  u.  JT 
CO  O  O 

I-  CO  -) 


js 

£  EJ  ©  O 
>  x  ©  C£ 
§  °  o  Q. 
£  -D  J2  fy 

co  a>  Q-fr; 

2  x 

I?  Oh 

.Q  ©  ©  ,_ 

co  o  o  O 
-J  O  O  2 

O  O  O  O 


©  E 

&  2 

2  * 

3  CO 

o  tr 

-C  CO 

<M  i o 

c  ° 
•=  o 

co  o> 

I  § 

o 

j:  §  -* 

|  JO 

£»  co  CM 

•  "O 

>*  o>  c  < 

•  8  "■ 
§!fi 

J  C  tt: 

**-  5  -o  ■ 
°  £  § 
8  15 

w  c 

2 1§  ■ 
.2  Q.  c 

o 

.2  o  o 
^  o  o  - 

jj  CO  CM  h 

(0004 


CO 

■o  _ 

©  m 

0)  g 

XI  k. 
^  3 
CO  o 

©  o 
o 

®  c 
2  .2 
5  c 

Q.  d) 

©  > 
CD  k- 

0) 

L.  4S 

1! 

o)  3 

c  © 

8  <D 

i  & 

J>*TJ 

©  2 


Q.  CD 

©  _c 
©  "■* 
j-  © 

a>  © 

1 1 
If 

£  O 
h-  2 

OO 


=*■0 

5S 

©  Q. 
©  © 
o  o 
o  o 

3  © 


E  © 

%  * 

>  © 

CO  tl 

■C  © 
©  CL 

*D 
*D  O 
©  O 
-Q  CD 

©  © 


..  O  O 

&  #  # 

-=  o 

2  00  a> 
©  o  o> 

(TOO 


Conclusion:  RUS  can  provide  stunning  improvements  i 
selected  areas  of  production  NDT 
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Initial  and  Final  measurements 
South  Girder 


2  2.5  3  3.5  4  4.5  5 
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We  have  not  succeeded  in  answering  all  of  our  questions. 
Indeed,  we  sometimes  feel  that  we  have  not  completely 
answered  any  of  them.  The  answers  we  have  found  only 
served  to  raise  a  whole  new  set  of  questions.  In  some 
ways  we  feel  that  we  are  as  confused  as  ever,  but  we  think 
we  are  now  confused  on  a  higher  level,  and  about  more 
important  things. 


-Author  unknown 
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Anderson  localization  °  Universal  conductance  fluctuations 

Ah  ronov- Bohw  effect  •  Mortal  elcclYM  persistent  currents. 


H?losi£E,c.  “  Phase  coherence  an  the  reelc  of  microns 


May  sWy  ^ihc- Ae portent  or  non-li near  effect 
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Penrose  “tile 


1791 


Bloch's  theorem  Statistics 
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Spectrum 


Density  of  states 
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PERIODIC  SUBSTRATE  —  1.1K,  7.5  layers 
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Simplified  Acoustic  Levitation 


Assume  only  two  forces  act  on  a  gas  bubble: 

1)  an  acoustic  force 

2)  buoyancy  force 


Also  assume  that  any  difficult  problems  can  be  ignored. 


Acoustic  Force  (Arb.  Units) 

At  equilibrium:  .1-0  .0.5  o.o  o.s 


acoustic^t  —  “  ^Fbuoyancy^t 
Facoustic  =  -V(z,t)VP(z,t) 

Fbuoyancy  =  Pl  V(z,t)  g  (pL  *  pL(t))  , 

4 

i 

P(z,t)  =  Pm  +  PAcos(kz)cos(cot)  : 

4 

1 

V(z,t)  =  Vm  -  VA(PA,co)cos(kz)cos(cot) 

t 

< 

Facoustic  =  -[Vm  -  VA(PA,CO)cOS(kz)COS(COt)]  < 

x[-kPAsin(kz)cos(cot)]  c 

<Facoustic>t  =  -(1/2)  kPAVAsin(kz)cos(kz) 

=  -(1/4)  kPAVAsin(2kz) 

Fbuoyancy  =  pL[V m  -  VA(PA,C0)cos(kz)cos(C0t)]g 
<Fbuoyancy>t  =  PL  Vm  g  (for  linear  oscillations) 
(1/4)  kPAVAsin(2kz)  =  pL  Vm  g 


Bubble  driven  below 
resonance. 


sin(2kz)  =  4  pL  Vm  g/(kPAVA)  =  4  pL  Vm  g/(kf(PA)) 
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Simplified  Bubble  Dynamics 


Rayleigh  -  Plesset  Equation 


Bubble  of  instantaneous  radius  R 
in  an  incompressible  fluid. 


Apply  Conservation  of  Energy 

KE  +  PE  =  constant 
KE  =  -  PE 


Kinetic  Energy 


KE  =  \'1phu1dV 

R 


Consider  the  mass  flux  through  a  spherical  surface  cenetered  on  the  bubble. 
If  the  fluid  is  incompressible,  then 

4  m2u{r)  =  4  7tR2u(R). 


So, 


u(r)  = 


u(R)R 2 
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Rayleigh  Plesset  Equation 


KE  =  j- 


Pl 


RR 


l2 


4nr2dr 


„  .2n4  7  dr 

=  2npLR  R  J  — 

Rr 

=  2npLR3R2 


R 


2 


Recall  that  the  effective  mass  of  a  small  (kR  «1)  bubble  is 


So, 


m«r=3(j“R3PL)- 


KE  =  \m,fR1. 


Taking  the  time  derivative  of  the  kinetic  energy  gives 

KE  =  2npL  [R3R2] 

=  2npL[2R1Ri  +2R3RR\ 


Rayleigh  -  Plesset  Equation 


Potential  Energy 


poo  =  pressure  in  the  water  far  from  the  bubble 
pL  =  pressure  on  the  “wet”  side  of  the  bubble  wall 


If  Poo  pL,  the  bubble  will  change  volume. 

The  time  rate  of  change  of  the  potential  energy  is 

PE  =  -(^)f 

=  ~{Pl  ~P„)4kR2R. 
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Rayleigh  -  Plesset  Equation 


Using 

KE  =  -PE 

gives  the  Rayleigh  -  Plesset  equation 


where 

2d  .  R 
Pl  =Pg+Pv-  —  ~aP-^ 


P  o  P acoustic 


Refinments  -  Keller  Equation 


Adding  the  effects  of  compressibility  leads  to  the  Keller  Equation 


The  SOI . nd  nclcl' 'i^cirivc,^ wt U,n ^1X7  ^1  "T- !°  8C"C''a,c  a,,(1  ohscrvc  SL. 
cmi.'cd  Ughl  iS  deeded  will,  a  VM  vZ^y  «he^I!Zl^  HIV) ' 


FIG.  18.  Simultaneous  plots  of  the  sound  field  (top),  bubble  radium  rm 
/=  22  3  kHz  m'neSC",C'  (b0"0m)  in  GLV21  P<  -  1*2  atm  a 


D.F.Oailan.  et.  a!.,  J.  Acoujl.  Soc.  Am.  91  ,  3  166-3  183  lig.m 
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B.  P.  Barber  and  S.  J.  Puttcrntan,  Nature  352.3ttt.32o  ,.< 


Occurrence  (%) 


Pulse  height  (V) 


0.5  0  54 


p  (atm) 


?5'  3  a;  P^lse  height  distribution.  A  pulse  height  of  0.1  V  represents 
11X10  photons  emitted.  20.193  kHt.  p.  Photons  per  burs,  as  a 
function  of  acoustic  pressure  amplitude.  /,  =  10.736  kHz 


B.  P.  Barber  and  S.  J.  Pulterm 


an.  Nature  352,318-320  (1991) 


Voltage  (mV) 


FIG.  2  The  average  of  single-pulse  outputs  of  the 
photomultiplier  tube.  a.  St  data  as  recorded  by  a 
conventional  (R928)  PMT;  b,  same  as  a.  but  using 
a  34-ps  laser  pulser  (Hamamatsu  PIP-01)  as  the 
light  source,  a  and  b  were  obtained  by  running 
the  PMT  output  into  a  digital  sampling  scope 
(HP  54 201  A),  c  and  d,  Data  for  a  microchannel- 
plate  PMT  (Hamamatsu  R1564U)  running  into  a 
20-GHz  digital  sampling  scope  (Tektronix  11802). 
c  is  for  the  St  source,  and  d  for  the  34-ps  laser 
pulser. 


B.  P.  Barber  and  S.  J.  Pullcrman,  Nature  352,318-320  (1991). 


♦  ■mo  (soc) 

PIG  I.  Voltage  versus  lime  at  (he  output  or  the  I’MT  Tor  SI.  nnd  the  IM.P. 

ic  imcse.i  cs  arc  chosen  so  that  the  two  curves  paSs through  the  50%  level 
a  the  same  lime  Ihcsc  curves  correspond  to  the  recording  or  about  25 
photoelectrons.  AHer  passing  through  the  delay  line  and  20-dI)  attenuator  a 
single  photoclcctron  corresponds  to  a  peak  amplitude  or  2  mV 


B.  P.  Barber,  ct.  a!.,  J.  Acoust.  Soc.  Am.  91  ,  306 1 -3063  (1992) 
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FIG.  3.  Histogram  of  events  versus  period  between  Hashes  lor  sonolumincs 
ccnce. 
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FIG *J  •  Ca,ibratcd  spectral  density  of  the  synchronous  pi¬ 
cosecond  flashes  of  sonolumincscencc  at  22  °C.  The  average 
spectral  energy  density  of  a  single  flash  can  be  obtained  by  di- 
vid.ng  by  the  acoustic  frequency  of  27  kHz.  The  dotted  line 
was  obtained  via  the  D  lamp  calibration.  The  points  with  error 
ars  were  obtained  by  calibrating  our  apparatus  with  a  QTH 

standard  of  spectral  irradiancc.  The  solid  line  is  a  25000  K 
blackbody  spectrum. 


R.  Hiller,  ct.  al.,  Phys.  Rev.  Lett.  69,  1182- 1  184 


(1992). 


222 


FIG.  4.  Raw  data  for  the  spectral  density  of  SL  at  10  and 
22  C.  The  peaks  have  been  chosen  so  that  the  curves  have 
equal  area.  These  curves  have  not  been  corrected  for  the  fiber 
grating,  or  PMT.  The  grating  is  blazed  at  300  nm. 


R.  Hiller,  el.  al.,  Phys.  Rev.  Lett.  69,1182-1184  (1992). 
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Fig.  5.  Time  dependence  of  SL  from  a  pure  N2 
bubble  in  water  (with  N2  dissolved  at  150 
mmHg):  (A)  low  drive,  (B)  high  drive.  The  SL 
intensity  has  been  normalized  to  the  emission  of 
an  air  bubble  at  the  standard  parameters  delin¬ 
eated  in  Fig.  1 .  Uncertainty  in  the  impurity  con¬ 
centration  is  about  0.05%.  The  long-term  mem¬ 
ory  (over  100,000  cycles  of  sound)  displayed  in 
this  data  is  indicative  of  an  as  yet-  unidentified 
physical  process  that  is  an  essential  aspect  of 
the  transition  to  SL.  We  were  unable  to  observe 
steady  SL  from  a  single  N2  bubble.  The  average 
radius  also  drifts  on  the  same  time  scales  in 
these  regimes.  Because  of  this  nonsteady  mo¬ 
tion  and  weak  emission,  we  were  unable  to  ob¬ 
tain  a  spectrum  of  a  N2  bubble. 
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FIGURE  5:  Values  of  ihe  intensity  of  sonoluminescence.  sound  field  level  P,.  maximum  bubble 
radius  Rm  and  ambient  radius  Ro  as  a  function  of  water  temperature  for  a  trapped  bubble  of  air. 

The  number  N  of  photons  per  burst  (with  wavelength  greater  than  200nm)  is  measured  in  each  case 
near  the  maximum  achievable  value.  The  bars  are  the  ranges  of  intensities  calculated  from  the  shock 
wave  theory  when  the  uncertainty  in  the  experimental  input  parameters  is  allowed  for. 
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Fig.  1.  Simultaneous  LC  scattering  at  scattering  angles  of  47° 
(dark  curve)  and  53°  (light  curve). 
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F^6-  and  pulsed  laser-scattering  graph,  showingthe  outpi 
of  the  PMT  s  used  to  monitor  SL  (upper  trace)  and  scatteri 
(lower  trace)  as  functions  of  time.  The  laser  pulse  starts  at  -( 
The  modulation  in  the  scattered  intensity  corresponds  to  L 
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Comparison  of  Multibubble  and  Single-Bubble  Sonoluminescence  Spectra 
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Comparisons  of  the  spectral  characteristics  of  sonoluminescence  from  cavitation  in  bubble  fields 
(MBSL)  versus  cavitation  of  single  bubbles  (SBSL)  have  been  made  for  aqueous  solutions  under 
similar  experimental  conditions.  In  particular,  alkali  metal  chloride  solutions  exhibit  sonoluminescence 
emission  from  excited  state  Na  or  K  atoms  in  MBSL,  while  SBSL  exhibits  no  such  emission.  Since 
the  metal  ions  are  not  volatile,  participation  of  the  initially  liquid  phase  must  occur  in  MBSL.  Surface 
wave  and  microjet  formation  in  cavitaling  bubble  fields  versus  the  high  spherical  symmetry  ot  collapse 
of  an  isolated  bubble  may  account  for  the  observed  differences. 


PACS  numbers:  78.60.Mq,  43.25.+y,  47.40.Nm 

It  has  long  been  known  that  under  certain  conditions 
acoustic  irradiation  of  a  liquid  can  result  in  light  emis¬ 
sion,  a  phenomenon  called  sonoluminescence  (SL)  [1,2]. 
The  process  typically  involves  the  application  of  high  in¬ 
tensity  ultrasound  to  a  liquid  by  an  immersed  acoustic 
horn  driven  with  a  piezoelectric  transducer.  The  resulting 
cavitation-bubble  field  is  made  up  of  a  complex  distribu¬ 
tion  of  gas  and  vapor-filled  bubbles  of  various  equilibrium 
sizes  that  pulsate  at  various  phases  relative  to  the  driving 
acoustic  pressure  field.  The  bubble  dynamics  is  further 
complicated  by  interactions  with  neighboring  bubbles  [3] 
as  well  as  with  the  vessel  walls.  Depending  on  the  loca¬ 
tion  within  the  pressure  field  and  these  other  influences, 
some  of  the  bubbles  may  grow  dramatically  during  the 
negative  portion  of  the  sound  field,  followed  by  a  quasi- 
adiabatic  collapse  that  results  in  the  heating  of  the  bubble 
interior  and  the  subsequent  emission  of  light  [4]. 

In  spite  of  the  complexity  of  cavitating  bubble  fields, 
many  studies  have  been  made  of  multibuljble  sonolumi¬ 
nescence  (MBSL)  and  the  influences  of  fluid  and  gas 
properties.  The  optical  spectra  of  MBSL  typically  con¬ 
tains  distinct,  pressure  broadened  molecular  or  atomic 
emission  bands.  Of  particular  significance  here  is  the 
identification  of  individual  transitions  from  excited  slates 
of  diatomic  carbon  (C2)  that  contribute  to  the  optical  spec¬ 
trum  of  MBSL  in  nonaqueous  liquids.  The  fitting  of  the 
measured  spectrum  of  C2  permitted  the  measurement  of 
an  effective  rotational  and  vibrational  temperature  of  the 
excited  states  of  C2  of  5100  K  [5]. 

Recent  experimental  advances  [6]  have  also  made  it  pos¬ 
sible  to  examine  both  the  temporal  and  spectral  nature  of 
sonoluminescence  from  a  single  bubble  (SBSL).  Here  a 
single  bubble  is  acoustically  levitated  in  an  aqueous  so¬ 
lution  that  has  been  partially  degassed.  The  bubble  can 
be  made  to  undergo  large-amplitude,  nonlinear,  presum¬ 
ably  radial  pulsations  during  which  light  emission  can  oc¬ 
cur.  Some  properties  of  SBSL  include  [7]  the  synchronous 
emission  of  light  with  each  and  every  acoustic  cycle,  tem¬ 


poral  flash  widths  of  less  than  50  ps,  and  a  continuous  spec¬ 
tral  energy  density  that  increases  lrom  the  visible  into  the 
UV,  with  eventual  fail  off  due  to  UV  absorption  by  the 
surrounding  water.  In  addition,  unlike  in  MBSL,  there  are 
little  or  no  electronic  or  molecular  bands  associated  with 
SBSL  spectra.  The  shape  of  the  spectrum  of  SBSL  has  lec 
some  researchers  to  suggest  that  SBSL  is  much  “hotter 
than  MBSL,  reaching  temperatures  as  high  as  50000  K 
[8],  and  possibly  much  higher  [9,10]. 

In  order  to  probe  the  differences  between  MBSL  am 
SBSL,  we  have  explored  emission  from  identical  aqueou: 
solutions  containing  potentially  emissive,  but  nonvolatile 
solutes  using  the  same  spectrometer  lor  both  systems 
Nonvolatile  solutes  provide  a  test  of  the  involvement  o 
the  initially  liquid  phase  surrounding  the  cavitating  bubbl 
in  the  sonoluminescent  event  [1 1,12].  An  observation  o 
an  SL  emission  peak  from  a  nonvolatile  solute  require 
either  that  a  fluid  shell  surrounding  the  bubble  be  heate 
sufficiently  [13],  or  that  liquid  droplets  containing  th 
nonvolatile  species  become  entrained  and  heated  withi 
the  bubble  [4,14]. 

It  was  not  possible  to  generate  both  MBSL  and  SBS 
in  a  single  apparatus,  owing  to  the  differences  in  tecl 
niques  used  to  generate  cavitation-bubble  fields  an 
isolated  single  bubbles.  MBSL  is  generated  using  a 
ultrasonic  horn,  which  produces  large  peak  pressure 
within  the  liquid  (around  10  bars),  with  high  levels  ( 
gas  saturation;  while  SBSL  occurs  with  applied  acoust 
pressure  amplitudes  near  1  bar,  and  gas  concentratic 
levels  that  are  a  fraction  of  saturation.  Nonetheles 
comparison  of  the  optical  spectra  can  be  made  by  usii 
identical  fluid  preparation  schemes  under  identical  gas* 
with  a  single  calibrated  spectrometer. 

The  SBSL  apparatus  consisted  of  a  quartz  cylindric 
levitation  cell  (8  cm  tall  by  4.5  cm  diameter),  as  shown 
Fig.  1  (a).  The  cell  was  closed  on  lop  with  a  glass  plate, 
hollow  cylindrical  PZT  transducer,  cemented  to  the  gku 
was  used  to  drive  the  levitator  in  tandem  with  a  pow 
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FIG.  4.  Chaotic  behavior  in  flash  variation  A I.  (a)  is  the  at¬ 
tractor  reconstructed  from  the  flash  lime  series,  (b)  is  the  dis¬ 
tribution  of  pointwise  dimensions  for  (a),  (c)  shows  the  results 
of  calculation  of  the  Lyapunov  exponents  for  the  attractor  in 
(a). 
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Observation  of  a  New  Phase  of  Sonoluminescence  at  Low  Partial  Pressures 
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The  acoustically  driven  pulsations  of  a  gas  bubble  lead  to  1 0 6-fold  changes  in  its  volume  and  the 
emission  of  a  light  flash  upon  collapse.  Mass  diffusion  between  the  bubble  and  the  gas  dissolved  in  the 
surrounding  fluid  maintains  this  steady-state  bubble  motion  only  at  low  partial  pressures,  around  3  Torr. 
This  diffusion-controlled  regime  is  uniquely  favorable  to  sonoluminescence  (SL)  from  hydrogenic  gases 
and  polyatomic  gases  with  low  adiabatic  heating.  Our  analysis  indicates  that  the  previously  investigated 
SL  from  bubbles  at  200  Torr  requires  a  nondiffusive  mass  flow  mechanism. 

PACS  numbers:  78.60.Mq 


A  gas  bubble  trapped  in  water  can  transduce  the  energy 
of  a  macroscopic  sound  field  down  to  the  microscopic 
level  where  it  is  emitted  as  picosecond  flashes  [1]  of 
ultraviolet  light  [2].  This  phenomenon,  sonoluminescence 
(or  SL),  is  particularly  sensitive  to  the  amplitude  of  the 
imposed  sound  field,  the  ambient  temperature  of  the  water 
[3],  and  the  gas  composition  of  the  bubble  [4].  This 
last  parameter  depends  on  the  concentration  of  the  gas 
dissolved  in  the  water.  In  particular,  the  observation  of 
SL  from  a  single  bubble  of  air  requires  that  the  water 
be  somewhat  degassed  [5],  whereas  with  pure  argon  the 
intensity  of  the  light  emission  is  relatively  independent  of 
the  dissolved  concentration  [4], 

Here  we  report  the  observation  of  a  new  phase  of  SL 
which  makes  this  phenomenon  accessible  to  hydrogenic 
and  polyatomic  gases.  As  shown  in  Fig.  1,  this  phase  is 
characterized  by  very  low  concentrations  of  dissolved  gas 
(about  10  ppb  for  deuterium)  or,  equivalently,  low  par¬ 
tial  pressures  of  solution  of  the  gas  in  the  water.  In  this 
region  of  parameter  space,  the  steady-state  motion  of  ev¬ 
ery  trapped  bubble  is  accompanied  by  light  emission,  as 
in  Fig.  2,  which  shows  the  stable  dynamics  of  an  ethane 
bubble  as  a  function  of  the  acoustic  drive  level.  Below  the 
lowest  drive  level,  the  bubble  is  unstable  against  dissolu¬ 
tion  in  the  water,  and  when  driven  at  an  amplitude  above 
the  upper  threshold  the  bubble  disappears.  We  have  also 
found  that  the  stability  of  light  emission  from  pure  noble 
gas  bubbles  dramatically  improves  as  the  partial  pressure 
is  reduced  to  the  level  where  light  emission  from  poly¬ 
atomic  and  hydrogenic  gases  is  optimized. 

The  investigation  in  this  region  of  parameter  space  was 
motivated  by  considering  the  mass  flow  between  the  pul¬ 
sating  bubble  and  the  gas  dissolved  in  the  surrounding 
fluid  [6-8].  When  the  bubble  expands  in  response  to 
the  rarefaction  of  the  driving  pressure,  the  pressure  of  the 
gas  inside  it  decreases  and  mass  diffuses  into  it  from  the 
fluid;  when  the  bubble  collapses,  the  gas  pressure  inside  it 
causes  mass  to  diffuse  out.  Requiring  the  net  mass  flow 
per  acoustic  cycle  to  vanish  specifies  the  steady-state  solu¬ 
tion  of  mass  diffusion  and  the  bubble  dynamics.  For  suf¬ 
ficiently  nonlinear  bubble  motion,  diffusive  equilibrium 


requires  that  the  partial  pressure  of  gas  dissolved  in  the 
fluid  be  given  by  [7,8] 

Pco/Po  -  3 (1) 

where  Px  is  the  partial  pressure,  P0  is  the  ambient  pressure 
(1  atm),  R0  is  the  ambient  radius  where  the  pressure  of 
the  gas  inside  the  bubble  is  the  ambient  pressure,  and  Rm 
is  the  maximum  radius  to  which  the  bubble  expands  in 
response  to  the  drive.  The  relation  (1)  is  derived  from 
coupling  the  measured  radius- versus-time  curves,  R(t ),  for 
driven  bubble  to  the  diffusion  equation  for  gas  dissolved 


FIG.  1.  Intensity  of  SL  for  ethane  and  deuterium  as  a  function 
of  partial  pressure  of  the  gas  dissolved  in  the  water.  The 
intensities  are  normalized  to  air  at  150  Torr.  Air  and  xenon- 
doped  nitrogen  exhibit  broad  peaks  centered  near  200  mm. 
Inset:  The  intensities  of  D2  and  C2H6  peak  at  partial  pressures 
of  a  few  Torr.  Between  50  and  150  inm  of  partial  pressure, 
ethane  gives  intermittent  light  with  an  intensity  of  about  0.03 
for  a  time  scale  of  less  than  25  s.  These  data  points  correspond 
to  light-emitting  states  stable  for  over  1  min.  The  experiments 
were  carried  out  in  the  acrylic  resonators  ( ojq  ~  35  kHz) 
described  elsewhere  [4J.  The  saturated  molar  solubilities  of 
D2,  C2H6,  xenon,  and  nitrogen  in  water  are  about  15,  41,  88, 
and  13  ppm,  respectively  (at  1  atm).  The  two  runs  with  ethane 
give  a  measure  of  the  systematic  variations. 
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The  spectrum  of  sonoluminescence  emitted  by  single  bubbles  of  Hi,  D2,  Hc\  and  He4  trapped  in 
H20  and  D20  has  been  measured.  We  find  that  heavy  water  has  a  dramatic  effect  on  the  spectrum 
of  hydrogenic  gases,  yielding  a  blackbody-typc  spectrum  with  a  spectral  peak  at  about  400  nrn.  The 
explanation  of  why  such  a  small  change  in  the  driving  fluid  leads  to  such  a  large  spectral  shift  is 
unknown. 


PACS  numbers:  78.60.Mq,  43.35.+d 

The  transduction  of  sound  into  light  (sonoluminescence 
or  “SL”)  is  medicated  by  the  extremely  nonlinear  pulsa¬ 
tions  of  a  trapped  bubble  of  gas  in  water.  Measurements 
indicate  that  the  huge  energy  focusing  involved  in  this 
process  is  accompanied  by  a  number  of  unknowns.  These 
include  the  light  emitting  mechanism  [  I  — 4J,  the  origin  of 
the  upper  and  lower  acoustic  drive  levels  which  delineate 
the  SL  regime  [5],  the  mechanism  of  mass  exchange  be¬ 
tween  the  bubble  and  gas  dissolved  in  the  liquid  [6,7], 
the  short  duration  of  the  flashes  [8,9]  (especially  their  fast 
turnoff),  and  the  limits  of  energy  concentration  that  can 
be  achieved  (i.e.,  what  is  the  spectrum  beyond  the  cutoff 
of  water?).  Finally,  why  is  water  the  friendliest  liquid  for 
,  SL  [10]? 

;  The  acoustic  properties  of  the  host  liquid  and  the 
‘  gas  bubble  are  obviously  key  physical  parameters  for 
]  this  phenomenon.  Thus  to  learn  about  SL  we  have 
!  measured  the  effect  of  isotopic  substitution  [11]  for  the 
gas  and  for  the  host  liquid.  In  this  direction  we  have 
observed  SL  from  H2,  D2,  He3,  and  He4  bubbles  driven 
by  sound  in  H20  and  D20  host  liquids.  In  view  of  the 
(also  unexplained)  strong  sensitivity  of  SL  to  the  fluid 
temperature  [5,12]  these  experiments  were  carried  out 
near  room  temperature  and  the  freezing  point,  which  for 
D20  is  3.8  °C. 

Figure  1  shows  the  cold  spectra  of  the  stable  helium 
isotopes  in  light  and  heavy  water.  A  general  feature  of 
these  measurements  is  that  the  spectrum  of  the  light  emit¬ 
ted  by  a  bubble  of  He3  is  slightly  more  ultraviolet  than 
that  of  a  bubble  of  He4.  Another  feature  is  that  bubbles 
in  heavy  water  are  somewhat  dimmer  and  less  strongly 
peaked  in  the  ultraviolet  than  are  bubbles  in  light  water. 
This  effect  of  heavy  water  on  SL  is  particularly  dramatic 
for  hydrogenic  gases  as  shown  in  Fig.  2.  The  appearance 
of  a  cutoff  to  the  spectra  of  H2  and  D2  bubbles  driven  by 
heavy  water  is  pronounced  and  unexpected,  and  stands  in 
marked  contrast  to  previous  spectral  investigations  of  SL 
[12].  The  rolloff  from  the  spectral  peak  occurs  at  the  peak 
efficiency  of  the  photodetector  [Hamamatsu  photomulti- 
I  plier  (PMT)  2027]  and  monochromator  (Acton  Research 
1  275,  1200  lines/mm  grating  blazed  at  300  nm). 
a  These  experiments  were  earned  out  in  a  cylindrical 
ig  quartz-walled  flask  (Heracus  supracil)  with  stainless  steel 


end  caps.  Resonance  frequencies  ranged  from  39.3  to 
40.6  kHz  for  heavy  water  and  41.2  to  42.6  kHz  lor  light 
water.  The  water  was  degassed  of  air  and  the  gas  of 
interest  dissolved  in  to  a  partial  pressure  of  3  mm  Hg 
in  the  water.  The  experimental  procedure  |I3|  involves 
seeding  a  bubble  with  a  Nidi  route  boiler  [14]  into  the 
scaled  flask  while  driving  the  flask  at  its  acoustic  res¬ 
onance.  The  bubble  was  aligned  along  the  optic  axis 
of  the  monochromator  through  the  use  of  a  laser  beam 
trained  on  the  output  slit.  Filters  were  used  to  suppress 
second  order  diffracted  light.  The  monochromator  slits 
were  3  mm  for  greatest  throughput  and  result  in  a  spec¬ 
tral  resolution  of  about  10  nm  FWHM.  The  radiance 
data  has  been  corrected  for  the  PMT,  monochromator, 
and  the  sapphire  window  of  the  cold  box.  The  data  have 


Wavelength  (nm) 

FIG.  I.  Spectral  radiance  of  helium  isotopes  in  light  and 
heavy  water  near  the  freezing  point.  Partial  pressure  of  gas 
is  3  mm.  Previous  measurements  of  helium  spectra  (131  were 
carried  out  at  150  mm  partial  pressure. 
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FIG.  2.  Time  evolution  of  case  (I)  from  the  nonadiabatic 
calculation,  p  is  in  kgm~J,  t>  in  kms“\  p  in  Pa.  and  T  in  K. 
Shown  are  solutions  for  (a)  t -20.490  257  (b) 

lb  “/c  +  0.156  ns;  (c)  ir  -/.+0.I9864  ns;  (d)  ij +0.19898 
ns;  (e)  +0.21641  ns;  (f)  t/-t„  +  0.222 35  ns;  (g) 

it  “'.+0.261  72  ns;  (h)  r*-/„  +  0.377  ns;  (i)  /,  +  0.396  ns; 
and  (j)  //-/.+0.4I7  ns. 


Wit  and  Roberts.  Phys.  Rev.  Lett.  70,3424-3427  (1993). 
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Time 

M(i.  .1.  Case  (I):  nonadiabatic  solution.  Tn  (»).  flic  bubble 
rsuliiis  and  the  shock  locations  (dashed)  arc  plotted  as  functions 
of  lime:  the  times  —  ij  employed  in  Tig.  2  arc  marked.  In  (b), 
(c).  and  (d),  the  maximum  density,  the  temperature  at  max¬ 
imum  density,  and  the  luminosity  of  the  bubble  arc  plotted  as 
functions  of  lime. 


Wtt  and  Roberts,  Phys.  Rev.  Lett.  7Q ,  3424-3427  (1993). 
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Is  sonolurninescence  collision-induced  emission? 

Lothar  Frommhold 

Physics  Department ,  University  of  Texas  at  Austin ,  Texas  78712-1081 


Anthony  A.  Atchley 

Physics  Department ,  Naval  Postgraduate  School,  Monterey,  CA  93943 
(August  16,  1994) 

An  estimate  is  attempted  of  the  collision-induced  emissioa  (CIE)  intensity  and  spectral  profile 
in  the  visible  and  near  UV  region  of  the  spectrum  of  N,-X  pairs,  where  X  represents  another  Nj 
molecule  or  an  argon  atom,  etc.,  for  conditions  that  correspond  to  shock  waves  believed  to  exist 
m  sonolurninescence  experiments.  Calculated  profiles  consist  of  superimposed  high  overtone  bands 
and  resemble  measured  profiles  of  sonolurninescence  spectra.  The  intensities  calculated  on  the 
basis  of  a  few,  simple  assumptions  concerning  the  induced  dipole  surfabe  compare  favorably  with 
measurements.  The  agreement  obtained  suggests  that  collision-induced  emission  is  a  viable  or  even 
an  attractive  alternative  to  bremsstrahlung  to  explain  sonolurninescence.  According  to  the  theory 

presented,  the  CIE  source  is  optically  thin  so  that  the  spectral  emission  profile  is  not  related  to 
Planck  s  radiation  law. 

43.20.+y,78.60.Mq,34.10.-fx,33.70.-w 


I.  INTRODUCTION 

A  single  bubble  of  air  may  be  trapped  in  an  acoustic 
standing  wave  setup  in  a  water  filled  container.  If  the 
proper  drive  frequency  and  amplitude  are  applied,  the 
bubble  may  emit  short  bursts  of  light,  an  effect  called 
sonolurninescence.  Some  recent  explanations  of  sonolu- 
minescence  have  centered  on  the  idea  that  a  converg¬ 
ing,  spherically  symmetric  shock  wave  is  launched  in  the 
gas  in  the  interior  of  the  bubble,  as  the  bubble  collapses 
under  the  influence  of  the  acoustic  standing  wave.  The 
assumption  has  been  that  the  temperature  in  the  (re¬ 
flected)  shock  wave  is  so  high  that  ionization  and  emis¬ 
sion  by  electronic  excitation  and/or  bremsstrahlung  oc¬ 
cur  (1).  Temperatures  in  excess  of  10*  K  are  often  as- 
sinned  in  such  work. 

Collision- induced  absorption  (CIA)  and  emission 
(CIE)  arise  from  dipoles  induced  by  intermolecular  in¬ 
teractions  (  collisions”).  Especially  the  absorption  spec¬ 
tra  of  the  common,  non-polar  gases  (e.g.,  hydrogen)  have 
been  studied  in  great  detail  in  a  number  of  laboratories; 
a  recently  published  monograph  summarizes  the  present 
knowledge  [2J.  Collision-induced  emission,  on  the  other 
hand,  is  less  well  studied  but  is  thought  to  be  an  im¬ 
portant  source  of  electromagnetic  radiation  in  the  atmo¬ 
spheres  of  planets  and  cool  stars  (3,4).  Infrared  emis¬ 
sion  of  shockwaves  has  also  been  explained  in  terms  of 
CIE  (5).  In  this  letter,  we  want  to  explore  the  possible 
connections  of  sonoluminenscence  and  collision-induced 
emission.  CIE  in  the  visible  from  shock  waves  at  tem¬ 
peratures  much  lower  than  10*  K  is  expected  at  high  gas 
densities. 

Caledonia  et  al.  reported  studies  of  interaction-induced 
light  emission  in  the  II2  fundamental  band  (near  2.4  /<m) 
at  densities  from  10  to  50  amagats  behind  reflected  shock 
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Measurements  of  collision-induced  emission  in  the  fundamental  vibration-rotation  band  of  hydro-, 
gen  are  presented  for  argon,  xenon,  and  neon  collision  partners.  These  absolute,  spectrally  resolved 
.  infrared  measurements  were  performed  at  high  densities  behind  reflected  shock  waves  over  the  tem¬ 

perature  range  of  900-3400  K.  The  emission  was  found  to  be  dominated  by  Q-branch  transitions  at 
high  temperature  due  primarily  to  the  dipole  moment  induced  by  the  overlap  between  the  electron 
clouds  of  the  collision  pair.  The  strength  of  this  interaction  was  evaluated  from  the  data  and  com¬ 
pared  with  similar  evaluations  determined  from  low-temperature  absorption  studies. 


I.  INTRODUCTION 

This  paper  reports  studies  of  the  coilision-induced 
emission  (CIE)  in  the  fundamental  vibration-rotation 
band  of  tlie  hydrogen  molecule.  The  spectrally  resolved 
infrared  emission  of  the  II2  band  was  measured  at  elevat¬ 
ed  temperatures  and  pressures  using  a  high-density, 
moderate-temperature  shock  tube  as  the  controlled  light 
source.  The  measurements  were  performed  for  several 
rare-gas  collision  partners  and  were  interpreted  within 
the  framework  of  existing  theory. 

Collision  induction  in  hydrogen  has  been  well  charac¬ 
terized  in  absorption  at  temperatures  below  400  K,1,2  but 
has  only  recently3,4  been  studied  at  higher  temperatures 
in  emission.  Collision-induced  processes  give  rise  to  vi¬ 
brational,  rotational,  and  translational  absorption  and 
emission  in  molecules  which  by  virtue  of  their  symmetry 
do  not  have  an  electric  dipole  moment  in  their  electronic 
ground  state.  The  proximity  of  a  collision  partner  makes 
possible  the  induction  of  a  transient  dipole  in  the  collid¬ 
ing  pair  resulting  from  interactions  due  predominantly  to 
the  permanent  quadrupole  moment  of  II2  and  to  the 
overlap  of  the  molecular  electron  clouds.  Such  interac¬ 
tions  produce  relatively  weak  radiation  signatures  which, 
however,  can  be  significant  at  high  pressures.  Indeed, 
CIE  by  H2  has  been  found  to  be  an  important  source  of 
infrared  radiation  in  planetary  atmospheres,5"7  and  its 
inverse  process,  collision-induced  absorption,  is  the  dom-* 
inant  source  of  infrared  opacity  in  stellar  atmospheres.* 

The  present  study  focused  on  the  high-temperature 
band  strength  and  spectral  shape  of  the  II2  fundamental 
vibration-rotation  band  centered  at  2.4  fi m.  The  band 
emission  was  measured  behind  reflected  shocks  over  the 
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temperature  range  of  900-3400  K.  Gas  mixtures  of  H2 
with  Ne,  Ar,  and  Xe  were  investigated.  Since  the  emis¬ 
sion  is  collision  induced,  its  intensity  scales  as  the  square 
of  the  density;  measurements  were  performed  over  the 
density  range  of  10-50  amagats. 

A  description  of  the  experimental  apparatus  and  tech¬ 
niques  is  provided  in  Sec.  II.  A  brief  review  of  the  theory 
of  Van  Kranendonk,9,10  and  how  it  was  applied  to  the 
present  data,  is  presented  in  Sec.  III.  The  data  are  sum¬ 
marized  in  Sec.  IV  and  analyzed  to  evaluate  the  induced 
dipole  parameters.  The  summary  and  conclusions  of  the 
study  are  given  in  Sec.  V. 

II.  EXPERIMENTAL  APPARATUS 
AND  PROCEDURES 

The  measurements  were  performed  in  the  Physical  Sci¬ 
ences  Inc.  (PSI)  shock  tube  facility.  The  shock  tube  is  a 
valuable  tool  for  this  type  of  study  because  it  can  readily 
heat  and  equilibrate  a  large  volume  of  gas  at  steady 
chemical  and  thermal  conditions.  Furthermore,  the  ther¬ 
modynamic  properties  of  the  gas  can  be  calculated  accu¬ 
rately  from  measurements  of  the  initial  gas  pressure  and 
shock  velocity.  On  the  negative  side,  a  measurement 
must  be  performed  during  a  relatively  short  test  time 
(  £  1.5  ms)  and  thus  there  is  little  opportunity  for  time 
averaging  the  signal. 

The  emission  measurements  were  performed  behind 
reflected  shocks  in  a  polished  stainless  steel  shock  tube 
with  an  interior  diameter  of  in.  The  driver  and  test 
sections  were  9  and  18  ft  long,  respectively,  with  |-in.- 
thick  walls.  A  schematic  diagram  of  the  apparatus  is 
shown  in  Fig.  1.  The  optical  measurements  were  made 
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FIG.  I.  Schematic  of  shock  tube  and  ir  instrumentation. 
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FIG.  1.  Total  emission  in  Watt  per  unit  source  volume, 
over  gas  densities  of  N*  and  X,  per  1  nm  wavelength,  for  eleven 
temperatures  from  6,000  to  100,000  K;  the  numbers  8  . . .  50 
are  short  for  8,000  IC  ...  50,000  K;  continuous  wave  emission 
is  assumed.  Also  shown  is  a  measurement  [18]  (circles).  The 
structures  discernible  at  low  temperature  are  artifacts. 


FIG.  2.  The  calculated  absorption  coefTicient  or  over  gas 
densities  of  N?  and  X,  for  eleven  temperatures  from  6,000 
to  100,000  K;  the  numbers  8,  ...  50  are  short  for  8,000  K 
...  50,000  K.  Structures  discernible  at  low  temperatures  are 
artifacts. 


Delay  Time  (ns) 


Average  of  1000  waveforms,  background  subtracted. 
Some  signal,  even  9  ns  after  main  bang. 

Jitter  in  system  is  =  5  ns,  but  not  yet  quantitative 
Preliminary  data  -  Is  there  a  preglow  as  well? 
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Sonoluminescence  as  Quantum  Vacuum  Radiation 
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Sonoluminescence  is  explained  in  terms  of  quantum  vacuum  radiation  by  moving  interfaces  between 
media  of  different  polarizability.  It  can  be  considered  as  a  dynamic  Casimir  effect,  in  the  sense  that  it 
is  a  consequence  of  the  imbalance  of  the  zero-point  fluctuations  of  the  electromagnetic  field  during  the 
noninertial  motion  of  a  boundary.  The  transition  amplitude  from  the  vacuum  into  a  two-photon  state 
is  calculated  in  a  Hamiltonian  formalism  and  turns  out  to  be  governed  by  the  transition  matrix  element 
of  the  radiation  pressure.  Expressions  for  the  spectral  density  and  the  total  radiated  energy  are  given. 
[S003 1-9007(96)00240-2] 

PACS  numbers:  78.60.Mq,  03.70.+k,  42.50.Lc 


Sonoluminescence  is  a  phenomenon  that  has  so  far  re¬ 
sisted  all  attempts  of  explanation.  A  short  and  intense 
flash  of  light  is  observed  when  ultrasound-driven  air  or 
other  gas  bubbles  in  water  collapse.  This  process  has 
been  known  for  more  than  60  years  to  occur  randomly 
when  degassed  water  is  irradiated  with  ultrasound  [I], 
Recently  interest  has  been  revived  by  the  contriving  of 
stable  sonoluminescence  [2,3]  where  a  bubble  is  trapped 
at  the  pressure  antinode  of  a  standing  sound  wave  in  a 
spherical  or  cylindrical  container  and  collapses  and  re¬ 
expands  with  the  periodicity  of  the  sound.  With  a  clock¬ 
like  precision  a  light  pulse  is  emitted  during  every  cycle 
of  the  sound  wave;  the  jitter  in  the  sequence  of  pulses 
is  almost  unmeasurably  small.  Shining  laser  light  upon 
the  bubble  and  analyzing  the  scattered  light  on  the  ba¬ 
sis  of  the  Mie  theory  of  scattering  from  spherical  ob¬ 
stacles  one  has  been  able  to  record  the  time  dependence 
of  the  bubble  radius  [4];  these  experiments  showed  that 
the  flash  of  light  is  emitted  shortly  after  the  bubble  has 
collapsed,  i.e.,  shortly  after  it  has  reached  its  minimum 
radius.  This  and  the  fact  that  the  spectrum  of  the  emit¬ 
ted  light  resembles  radiation  from  a  black  body  at  sev¬ 
eral  tens  of  thousands  degree  kelvin  have  led  to  the 
conjecture  that  the  light  could  be  thermal  radiation  from 
the  highly  compressed  and  heated  gas  contents  of  the 
bubble  after  the  collapse  [5].  It  has  also  been  argued 
that  the  experimentally  observed  spectrum  would  equally 
well  be  compatible  with  the  idea  of  a  plasma  form¬ 
ing  at  the  bubble  center  after  the  collapse  and  radiating 
by  means  of  bremsstrahlung  [6].  An  alternative  sugges¬ 
tion  has  tried  to  explain  the  sonoluminescence  spectrum 
as  pressure-broadened  vibration-rotation  lines  [7],  but  al¬ 
though  this  theory  has  been  very  successful  in  the  case  of 
randomly  excited  (multibubble)  sonoluminescence  seen  in 
silicone  oil  it  has  been  inefficacious  for  sonoluminescence 
in  water. 

All  of  the  above  theories  have  serious  flaws.  Both 
blackbody  radiation  and  bremsstrahlung  would  make  a 
substantial  part  of  the  radiated  energy  appear  below 
200  nm  where  the  surrounding  water  would  absorb  it.  If 
one  estimates  the  total  amount  of  energy  to  be  absorbed 


corresponding  to  the  observed  number  of  photons  above 
the  absorption  edge,  one  quickly  convinces  oneself  that 
this  would  be  far  too  much  to  leave  no  macroscopically 
discernible  traces  in  the  water,  as,  for  instance,  dissocia¬ 
tion  [8];  however,  nothing  the  like  is  observed.  Another 
very  strong  argument  against  all  three  of  the  above  the¬ 
ories  is  that  the  processes  involved  in  each  of  them  are 
far  too  slow  to  yield  pulse  lengths  of  10  ps  or  less,  but 
which  are  observed.  Moreover,  if  a  plasma  were  formed 
in  the  bubble,  one  should  see  at  least  remnants  of  slow 
recombination  radiation  from  the  plasma  when  the  bubble 
reexpands.  As  to  the  theory  involving  vibration-rotation 
excitations,  the  line  broadening  required  to  model  the  ob¬ 
served  spectrum  seems  rather  unrealistic. 

In  its  concept  the  theory  to  be  presented  here  has  been 
loosely  inspired  by  Schwinger's  idea  [9]  that  sonolumi¬ 
nescence  might  be  akin  to  the  Casimir  effect,  in  the  sense 
that  the  zero-point  fluctuations  of  the  electromagnetic  field 
might  lie  at  the  origin  of  the  observed  radiation.  More 
closely  related  to  this  is  the  Unruh  effect  well  known 
in  field  theory  [10];  its  original  statement  is  that  a  uni¬ 
formly  accelerated  mirror  in  vacuum  emits  photons  with 
the  spectral  distribution  of  blackbody  radiation.  However, 
the  phenomenon  is  far  more  general  than  that  and  in  par¬ 
ticular  not  restricted  to  perfect  mirrors.  This  kind  of  quan¬ 
tum  vacuum  radiation  has  been  shown  to  be  generated 
also  by  moving  dielectrics  [11].  Whenever  an  interface 
between  two  dielectrics  or  a  dielectric  and  the  vacuum 
moves  noninertially  photons  are  created.  In  practice  this 
effect  is  very  feeble,  so  that  it  has  up  to  now  been  far 
beyond  any  experimental  verification.  Sonoluminescence 
might  be  the  first  identifiable  manifestation  of  quantum 
vacuum  radiation. 

The  mechanism  by  which  radiation  from  moving  di¬ 
electrics  and  mirrors  in  vacuum  is  created  is  understood 
most  easily  by  picturing  the  medium  as  an  assembly  of 
dipoles.  The  zero-point  fluctuations  of  the  electromag¬ 
netic  field  induce  these  dipoles  and  orient  and  excite  them. 
However,  as  long  as  the  dielectric  stays  stationary  or  uni¬ 
formly  moving  such  excitations  remain  virtual;  real  pho¬ 
tons  are  created  only  when  the  dielectric  or  mirror  moves 
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WHAT  IS  NONLINEAR  ACOUSTICS? 

Sources  of  nonlinearity  in  fluids: 

1)  Equation  of  state 

2)  Convection 


Effects  on  plane  waves: 

•  Propagation  speed  varies  along  waveform 


Consequences: 

•  Wave  distortion 

•  New  frequencies 

•  Radiation  pressures 

•  Acoustical  streaming 

•  Shockwaves 
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Except 


Substance 
distilled  water 


Pressure 
1  atm  y 
200  kg/ cm 
4000 
8000 

sea  water 
(3.52) 
methanol 
ethanol 


n-propanol 

N-butanol 

acetone 

beneze 

chlorobenzene 
liquid  nitrogen 
benzyl  alcohol 
diethylamine 
ethylene  glycol 
ethyl  formate 
heptane 
hexane 


R.T  er 


TABLE  I 


Values  of  B/A. 

Where  Indicated,  All  Values  are  at  Atmospheric  Pressure 


ix  ,lq 

B/A 

Substance 

T,  °C 

B/A 

$ 

0 

4.2 

methyl  acetate 

30 

9.7 

20 

5.0 

cyclohexane 

30 

10.1 

40 

5.4 

nitrobenzene 

30 

9.9 

60 

5.7 

mercury 

30 

7*8 

80 

6.1 

sodium 

110 

2.7 

100 

6.1 

potassium 

100 

2.9 

tin 

240 

4.4 

30 

5.2 

indium 

160 

4.6 

30 

6.2 

bismuth 

318 

7.1 

30 

6.2 

30 

5.9 

monatomic  gas 

20 

0.67 

diatomic  gas 

20 

0.40 

20 

5.25 

methyl  iodide 

30 

8.2 

20 

9.6 

sulfur 

121 

9.5 

0 

10.4 

glycerol  (42  H_0) 

30 

9.0 

20 

10.5 

1,2  -  dichloro- 

30 

11.8 

40 

10.6 

hexafluoro- 

20 

10.7 

cydopentene 

20 

10.7 

(DHCP) 

20 

9.2 

20 

9.0 

30 

9.3 

b.p. 

6.6 

30 

10.2 

30 

10.3 

30 

9.7 

30 

9.8 

30 

10.0 

30 

9.9 

j  AJgwfi  n-gg-r 


Characterization  of  biological  madia:  L  Bj$m+ 

Table  1  Some  8/ A  values  for  biological  material* 


Biological  material  (and  state) 

Method* 

8/A  (and  uncertainty) 

Ref* 

1. 

Bovine  serum  albumin  (BSA)  (20  g/100  cm3,  25*C) 

Therm. 

6.23  (±  0.25) 

31 

BSA  [22  g/100  cm3.  30‘C) 

FA. 

6.45  (±  0.30) 

21 

BSA  (38.9  g/100  cm3.  30‘C) 

FA. 

6.64 

30 

BSA  (38  9  g/100  cm3,  30*C) 

Therm. 

6.68  (±  0.2) 

30 

2. 

Haemoglobin  (50%.  30*C) 

F.A. 

7.6 

22 

3. 

Whole  porcine  blood  (12%  haemoglobin,  7%  plasm* 
proteins.  30*C) 

FA. 

6.2  (±0.25) 

22 

4. 

Beef  liver  (Whole.  23‘C) 

F.A. 

7.75  (±  0.4) 

22 

Beef  liver  (Homogenized.  23‘C) 

F.A. 

6.8  (±0.4) 

22 

Beef  liver  (Whole.  30‘C) 

F.A. 

6.42 

30 

Beef  liver  (Whole.  30‘C) 

Therm. 

6.88 

30 

Beef  liver  (Whole.  30*C) 

Therm. 

6.54  (±  0.2) 

32 

Dog  liver  (30*C) 

F.A. 

7.6  -7.9  (±0.8) 

35 

Pig  liver  (25*C) 

FA 

6.7  (±1.5) 

36 

Human  Irver  (Normal.  30‘C) 

FA 

7.6  (±0.8) 

35 

Human  liver  (Congested,  30‘C) 

FA 

7.2  (±0.7) 

35 

5. 

Pig  fat 

Therm. 

10.9 

27 

Pig  fat 

FA 

11.0-11.3 

27 

Human  breast  fat  (22‘C) 

Therm. 

9.21 

32 

Human  breast  fat  (30‘C) 

Therm. 

9.91 

32 

Human  breast  fat  (37‘C) 

Therm. 

9.63 

32 

8. 

Canine  spleen 

FA 

6.8 

37 

Dog  spleen 

FA 

6.8  (±0.7) 

35 

Human  spleen  (Congested) 

FA 

7.8 

37 

Human  spleen  (Normal.  30‘C) 

FA 

7.8  (±0.8) 

35 

7. 

Beef  brain  (30‘C) 

FA 

7.6 

27 

8. 

Beef  heart  (30*C) 

FA 

6.8-74 

27 

9. 

Pig  muscle  (30‘C) 

FA 

7. 5-8,1 

27 

Pig  muscle  (25‘C) 

FA 

6.5  (±  1.5) 

36 

10. 

Oog  kidney  (Normal,  30*C) 

FA 

7.2  (±0.7) 

35 

Canine  kidney  (30*C) 

FA 

7.2 

37 

11. 

Human  multiple  myeloma  (22‘C) 

FA 

5.6 

32 

Human  multiple  myeloma  (30'C) 

FA 

5.8 

32 

Human  multiple  myeloma  (37#C) 

FA 

.  ** 

32 

•Therm.  ■  thermodynamic  method;  FA  *  finite  amplitude  method 
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the  stste  of  the  tissue.  Cancerous  tissue  normally  shows  a 
higher  water  fraction  than  normal  tissue.  The  water 
fraction  goes  from  0.76  in  normal  liver  tissue  to  0.90  (br 
multiple  myeloma”.  In  general  water  in  tissue  mav  be 
found  as  bound  water  and  as  free  water  in  equilibrium 
with  one  another  and  expressed  by 

(HaO)„  *  nH^O  (12) 

where  (H20)„  is  referred  to  as  bound  water  while  H,0  it 
generally  referred  to  as  free  water  An  increase  in  the 
bound  state  means  that  on  average,  molecules  have  a 
greater  degree  of  association  with  the  neighbouring 
molecules  which  means  that  they  are  held  more  strongly 
together.  This  stronger  binding  also  makes  a  larger 
ultrasonic  pressure  necessary  in  order  to  stretch  the 
intermolecular  bonds  into  their  non-linear  region,  which 
m  tcroscopically  is  being  felt  as  decreased  non-linearity  of 
the  water  according  to  Equation  (3).  This  suggests  that  the 
magnitude  of  B/A  in  water-like  media  may  be  related  to 
the  relative  amounts  of  bound  and  free  water. 

The  equilibrium  between  these  two  water  states,  for 
instance  expressed  through  the  ratio  of  bound  to  free 
water,  is  closely  related  to  the  state  and  the  nature  of  the 
tissue  as  shown  by  NMR  studies*  Prospective  relations 
between  the  ratio  of  bound  to  free  water  and  the  non¬ 
linear  parameter  B/A,  have  recently  been  suggested*  It 
was  concluded  by  Yoshizumi  et  at*9  that  the  temperature 
dependence  of  B/A .  of  water  for  instance,  could  be  due  to 


the  change  in  the  ratio  of  the  bound  to  the  free  waterwfc 
the  change  in  temperature.  Whether  the  estimates  of  tW 
ratios  of  bound  to  free  water  determined  from  B/A 
measurements  can  be  used  forcharacterizatioi  of 
biological  media  such  as  human  tissue  is  still  an  ooq  i 
Question  which  has  to  be  studied  more  closely.  ! 

The  existence  of  severalpossible relationships  betww*  I 
B/A  of  biological  media  and  other  physical  qualities  sack 
as  intermolecular  potentials,  macrostructure,  water  i 
fraction  and  ratio  of  bound  to  free  water  of  the  biological  I 
media  emphasize  the  need  for  further  systematic  studio.  ' 
This  will  probably  demand  an  internationally  funded  * 
research  programme  where  several  qualified  laboratories 
in  various  countries  share  a  research  programme  om 
several  years  on  advanced  modelling  of  biological 
media. 

In  spite  of  the  inhomogeneous  character  of  bioiogied 
materials  such  as  tissues,  where  scattering,  plus 
cancellations,  dispersion,  etc.  influence  the  ultrasonic  * 
wave  propagation,  a  reliable  experimental  precede* 
leading  to  reproducible  B/A  data  should  be  developed  » 
prove,  in  vitro,  that  B/A  may  be  used  for  the  charade 
ization  of  biological  media.  If  the  uncertainties  found  is 
relation  to  the  experimental  data  could  be  reduced  forte 
thermodynamic  or  for  the  finite  amplitude  method  -  <*  \ 
maybe  for  both  -  an  answer  to  the  question  of  appl** 
ability  of  B/A  for  characterization  of  biological  media  may  t 
be  found.  i 

The  prospective  development  of  a  clinically  applicable  ' 
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Figure  1.1 

Demonstration  o!  the  collineer 
modulation  of  sound  by  sound 
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G«m+j  2ni>,2>9  (11 71) 


(Nonlinear  Interaction  of  a  Ffnlte-Amolltude  Wave 
with  a  Small-Signal  Wave  In  Air) 


R«C.  SI 
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COEFFICIENT  OF  NONLINEARITY 
FOR  COMPRESSION  WAVES  IN  ELASTIC  SOLIDS 


Shock  Formation  Distance: 

x  =  p7k 

where 


e  =  Mach  number  or  strain 
k  =  wavenumber 
0  =  coefficient  of  nonlinearity 


/  3  A  +  3R +  C\ 

v2  +  ptf  J 


substance 

0 

air 

water 

iron 

plastic  foam 
marble 
sediment 
rock 

1.2 

3.5 

4.4 
~  100 
~  800 
~  102-103 
>  104 

L.  A.  Ostrovsky,  J.  Acoust.  Soc.  Am.  90,  3332  (1991) 
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Nonlinear  Energy  of  Deformation 


•  Isotropic  solid  to  cubic  order  in  strain  is  (Landau’s  notation) 

£  =  fi&ik  +  (2^  ~~  3^)  eii  2^c»*c*jc*i  d"  &eikejj  d-  ^ejj ’ 

where 

1  /  dut-  duk  duj  duj  \ 

**'*  =  2  +  ax' +  axj 

is  the  Lagrangian  strain  tensor.  Thus, 


E  = 
£2  = 

E3  = 
+ 


£"2  +  £3, 

1  /  dut-  du*  \ 

4F\dXk  dX, ) 


(  1  A  dui_duj_duj_  /  s  £  _  m  <t±  ( ita. 

V  +  4  /  dX,dXk  V2  2  3/  3X, 

3 u;  3ui  3uj  B  duj  duk  duj  C  /  &ui  V 
12  dXk  dX,  dX,  +  2dX~kdX,  dX,  +  3  \dXj 
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SHEAR  WAVES 


Linear  Surface  Waves 


•  Rayleigh  (1885) 

•  Stoneley  (1924) 

•  Scholte  (1942) 


Stoneley  Wave 


-x 

y 


Scholte 


II.  WEAK  SHOCK  THEORY 
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uieflK  sUe/tt  th&sRV 


v/iW=MelocK^ 

X 


Xi,(-t\  =  post -fro  A 


J  \j)gp(t  fkiotk  L'lMi'f  q£  flainlwe-  UuQQ<\lo4’,  ftcfaHoJIS. 


A/S=G+  ■^■(Ua.+  Ul,) 


•7T  /  ayjaa^  Equal  Area.  ftolg 


^£*tcc<\j\M£  pcsc-fVcift  X& 

e4  shecff  log 


areas 


(4i  =  $2- 


9X 
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sVcsU  s-fruc-fure  (c-<?..  rKe  4tiwe^  r\c+  described 
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CH.  1.  PLANE  WAVES  IN  MEDIA  WITHOUT  DISPEBIC* 


air 


FIg»  1*13.  Change  in  the  shape  of  a  single  tri¬ 
angular  pulse. 

Everything  that  has  been  said  in  this  section  applies  to  the 
behavior  of  a  periodic  signal  in  a  nonlinear  medium,  which  is 
harmonic  at  the  input  to  the  system.  As  a  conclusion,  we  consider 
tne  progress  of  die  propagation  of  a  uniform  disturbance. 

Let  the  profile  of  the  Initial  disturbance  have  the  form  of  an 
isoacelea  triangle  ABC,  drawn  in  Fig.  1.13.  Hie  discontinuity  is 
formed  here  at  the  point  <r  »  m0t/2.  S  Is  easy  to  see  that  in  the 
region  0  <  a  <  ojt0/2,  the  area  of  the  profile  does  not  change  ami 
is  equal  to  «r*/ 2.  When  the  poise  passes  the  point  o  -  wr»/2.  Its 

**“  po,““  “•  -*  >  "V*. 

As  we  already  know,  in  place  of  CB*  we  have  the  discontinuity 
ED,  constructed  so  that  the  areas  of  the  triangles  CEF  and  FDB* 
are  equal.  But  this  means  that  the  areas  of  the  triangles  ABC  aid 
ADE  are  also  equal,  i.e.,  the  conservation  of  momentum  holds  in 
lowest  approximation.  Ihe  formation  of  a  discontinuity  conse¬ 
quently  leads  to  the  spreading  out  of  the  disturbance  and  a  •smear¬ 
ing  out*  of  its  amplitude  throughout  the  medium. 

ut.  Kr^an  **•  <***«  of  equality  of  the  areas  ABC  aid  ADE, 

tan  °  = ,  we ean  eatabliah the  formulas  which 

describe  the  spreading  out  at  a  >  «To/2: 


(5-21) 


and  the  decrease  in  the  amplitude  of  the  discontinuity  ED  -  AE  tan  a: 

EDr-  1 


V+TT' 


which  aher 
periodic  v 
of  the  init 


(5.22) 
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N  WAVES 


Sources: 

•  sonic  booms 

•  sparks 

•  lithotripters 

•  thunder 

•  explosions 

•  bullets 


I 

Af  ~  >0T 

T  ~  vT/T 


(rise  time)  (duration) 
1 03  sec  1 0'1  sec 

1 0'6  sec  1 0'5  sec 

1 0-9  sec  1 06  sec 
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THEORETICAL  PREDICTIONS  OF  THE  ACOUSTIC  PRESSURE 
GENERATED  BY  A  SHOCK  WAVE  UTHOTRIPTER 

A.  J.  Coleman,  M.  J.  Choi  and  J.  E.  Saunders 

Medical  Physics  Department,  St.  Thomas’  Hospital,  London  SE1  TEH,  UK 
(Received  5  July  1990:  in  final  form  I  October  1990) 


Ultrasound  in  Med.  A  Biol.  VoL  17.  No.  3.  pp.  245-233.  1991 
Printed  in  the  U.S.A. 


Fig.  4.  Plots  of  the  peak  positive  pressure  gain  (P+/PJ  along 
the  beam  axis.  zlF,  calculated  for  a  pulsed  (exponentially 
damped  sinusoidal)  aperture  waveform  with  peak  pressures. 
p„  of  1.3, 5  and  7  MPa  (n^  a  384).  The  curve  for  5  MPa  is 
shown  as  a  dotted  line  and  corresponds  to  that  predicted  for 
the  Dornier  HM3  operated  at  around  20  kV 
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6curte  Condli'-tron- 

u(o}t)  =  U,  £aia  or 

Fourier  Series  Soiu-b’cnA  *. 

# 

U-CtTi  f)  =  Uo  ?  saw  rtw-v  o-  =• 

6a i  <r<  i  (fuU;f/X 

=  (sao+astl^ 


Fio.  4.  Variation  of  2nd  Harmonic  inaptitude  with  fre¬ 
quency  of  fundamental. 


Fic.  5.  Variation  of  2nd  harmonic  premure  with  funda 
mental  preeaure. 
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Use  s.q (a  solution  ter  ^.3) 
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FIGURE  1-1 
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FIGURE  4-1 

BLOCK  DIAGRAM  OF  THE  EXPERIMENTAL  APPARATUS 
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FIGURE  5-5 

AMPLITUDE  RESPONSE  CURVES  AT 
FUNDAMENTAL  FREQUENCY 


FIGURE  5-6 

AMPLITUDE  RESPONSE  CURVES  FOR 
THE  FIRST  THREE  HARMONICS  OF 
AN  INITIALLY  SINUSOIDAL  WAVE 
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FIGURE  3-8 

A  SCHEMATIC  DIAGRAM  OF  MODIFIED  WEAK -SHOCK  THEORY 
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EXPERIMENTAL 


COMPUTED 


Figure  l4  A  comparison  of  experimental  and  computed  time  waveforms 
for  a  wave  produced  by  an  800  Hz  source  at  one  end  of  a  96  ft 
progressive  wave  tube.  The  source  SPL  is  159  dB  (re  0.0002  pbar); 

VA  stands  for  the  vertical  amplification  of  the  oscilloscope.  The 
dashed-line  predictions  are  given  by  weak-shock  theory  alone,  that 
is  without  any  account  taken  of  tube  wall  attenuation  and  dispersion. 
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The  second  example  problem  waveform 
(a)  x  s  0,  (b)  0  <  x  <  x,  (c)  x=5 
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Figure  3-5 

The  second  example  problem  (cont  d) 

(o)  The  source  waveform  in  the  shifted  coordinate  system 
(b)  The  souree  waveform  after  propogoling  25  ft. 

(e)  Just  prior  to  shock  merger. 

(d)  The  simplified  waveform  after  shock  merger. 


FIGURE  6-15 

EXPERIMENTAL  SPECTRUM  AT  VARIOUS  DISTANCES 


D.  A.  Webster  and  D.  T.  Blackstock:  Collinear  interaction  of  noise  with  «  tone 


FIG.  8.  Computed  tone-noise  interaction  spectra.  Symbols: 
a&x/x,  x  =  7.45  m.  The  conditions  are  roughly  those  of  ex¬ 
periment  I. 


III.  MODEL  EQUATIONS 
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V/lSCgOju  HFftTeoNlbUCTfAlfr  FLU 
Cco-tiouHq '. 

^|  +  fV*u  =  0  (■fe’s=^  +  a*^) 


/Vc'<M€lrt,fU0 A  1 

^+^  =  L\*spZ  <*■  +(i+|-)^^ 


=  Fs.(f^f(§ir)  ,  perfet<ja<. 


fiee  Looaau  |  IrfsW.-ta ,  £!aii  M«haflic<A 
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HOb£L  £$UflT7gWi<;  Of  AtoMl/AlEflR  ACOUSTICS 

*  £xac£  ctjaa 4700  -for  b^sbsS  Cj^/K-o) 

per -feet  qas  in  -terns  bf  uelocrlw 

fc^-b'alCU:  3  J 

ci  ~4bt  ~  (2  +"2  ^ 

+  M(<£fc + i  lv<ttz)  v2^ 

*  Con^OiT*  S+d\r47VU»  rt-t-  i'O  a€rZ>e(as47C<4M  i 
anii  “fee  p-eHur(&o.4f<i<0  techniques  uiea  »u 
/loviline&r  aasus+t'cS. 

*  ftbMA/MTAGE  :  Xt’s  cKact 

*  bl S^bV 

/}  restricted  Id  /essl-ecs  *ases 

2^  no  ex-ad-  ^Mio^s  (-exc^px 

Pofssofl  So/u+7oW  -fer  pnxjfessivl'e  ploort^ 
uoviec^  u 


fHi  linear  -ferAis  Ce-4-«  V-p)  are 

OCt),  w^re  d 

1  it 

-  accasb'c  /Hactt  fluwta 


a 


A# 


•  frit  tas,  c©e'f£rc?evrfs  are  6(p) • 

(m) 

0  Rrsi-orcier  -tcrtMS  C 

6Lz)  j  (©&s(ess  ( i  near 

2\  S«<:oa^ -order  -terms'' 

OCl^j  lossless  ^uadra-h'c. 
0(jAi^  j  (o<,s,m  linear 

3)  Rfjh-er-ordeC  4rroA6*. 

6(t3)j  Q[ptz)  )  0 (pzt))  e-fcc. 

*••  bt'scard  alj_  ^ifg(r>er-ordier  iVl 

<den\roi4-icfl£  <sf  all  00C<i-e(  ^^ug-br'cnA  S 
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l 


u 


Let  f  =  P-Po  ,  f'=P-f0i  Use  ^rst-crde< 

c«!<xfioi»^  -fo  sryivtafrH}  rcr'WS  > 

^^4  lao^re  vcr-h'di^  Cvxu)  b  db+affl. 

r,ofl4:lAU>^ 

-^+f.v-u.-|rg3g:  + 


=  LacjraA<j?Q-'A 


Full  second- order  U'gue  e^u^-fifliA. 
[fortoAseo  etoiM  3£,  CH^yj 


Ul^s^enM-fc  eoua4-foft 

[U3esttrrvj€l-tj  JflS#  3£3  S3-S-  ((*1 4.3)3 


,  S  _ 

+  3  2-t*  ~  fe^  Z-t*- 


TKe  app©Y.iv\\a.+icrv»  -fo  o\o4tiI^  \j!tes>4ecud'fc 

eaua'Miru  revfcricf^  i-t  -fo  ^uas.ipl&'rt'fi 

vaJ a\ies.  1 


^  =  (  Hr  —  £c>e£&d€>fl+  of  Aov\(r<  f\ean4:H 
=  Sc  a  fid  <A'i'f^u^(\j'rfcy 

Cr*£- 

=  ^x  = -fKe-r«viovl  4ccu  s  a4±ev\ua4-fcru. 
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’  feurgers  £ao,oJri(n^ 

[KUUl(nJ-  et  fll. ,  {kusWfca  iff  ,  24  S 

BeGiA  uitl)  Uie«»-tervi€(t  e^ua-h'cn/i  -for  hb  ■ 

tefr  a- sty?  -  czJt1  $&W 

TLio  approxi'AkvU  sdCcs  -fee  liW-Hvid  cages’. 

-p^exp  )  £>  =  c> 

d  iif\(^/t  +  (^'SSP'fi) 

Bo4-U  &oluJ->O'0&  ■for'M 

-  _  -/--•$-  =  rc-bo.  oded.  -ti<v\£  7 

^(=  £)<  =  K<>te»\»i ^  (ev^bh  £ca-(€. 


a?  J1 

a-tf  -  zt-*- 
a*  l  2? 


a* 


a? 


o-  *  <->  -  g:  &  c~  - 

5><£  -  c|  Srt2-  CoVUh' t  T  Z.  *xf- 

({c-ep  OAl'j  S^C^S^d-OC^^r  £>C&k?' 


S£  4  A 

^  =i5^  'h^S5  at 
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spectral  moHefticfl-L  sqlqttoM 

t>i'(Na€rtS('cnAl«s,S  Supers  C^aa-hoiA: 

*p  _  i  ££ 


p*-2- 


+  p 


A  _ 

P=^  j  ^"=X  )  * 


i 

aR 


Foon'er  Scries  Ey^pairtSfcrvA 


PCr,+)=i  21  TU<»Wor 

a--N/ 

R«sa/*t('n^  Coopted  f  (Ysrt-Crd-er  CbE's  - 

A  f 

—  __  -D_  p  1^  + 

-  p  +  4-  Vf^i 


Pa/\  ^ 


(VA=<l-hj 
J  \ _ 


(NA  1  <V\ 


+) 


a.fciorpVid'A  siw*  ^ <■('€■  "ffc^s  • 

Smpte  <rv\  fey'arioi'bro-r^  aJoscrphW 

arxji  disPtrUcvi 


jv 

r 


Aa  -f-j  b 


316 


******** 


********  Choose  absorption  law  (ABSTV,  ABSREL  or  ABSBL)  for  A(N) 

DO  10  N=1,NHAR 

P(N)  =  CMPLX(0.,0.) 

A(N)  =  ABSTV (N) 

10  CONTINUE 

P(l)  -  CMP LX (0 . * -1 • ) 


********  Second  order  Runge-Kutta  ******** 


20 

30 


40 


00  KIW^STEP  *  (  (J*N/4. )  *  SUM(P,N,NHAR)  -  A(N)  *  P  (N)  ) 
P2(N)  =  P(N)  +  K1(N) 

CONTINUE 

00  K 2  =_ STEP^*  (  (J*N/4. )  *  SUM  (P2,N,  NHAR)  -  A(N)  *  P2  (N)  ) 
p  (N)  =  P(N)  +  (Kl(N)+K2)/2. 

CONTINUE 


X  =  X  +  STEP 

IF  (X.LT.XOUT (I) )  GO  TO  20 


DO  50  N=l, NHAR 

POUT (I, N)  -  P(N) 

50  CONTINUE 

1  =  1  +  1 

IF  (I.LE.IOUT)  GO  TO  20 

********  Write  complex  harmonic  amplitudes  to  file  "Harmonics"  ******** 


DO  60  N=1,NHAR  ,  T , 

WRITE (12,*)  N,  (  POUT (I, N),  I-l,IOUT  ) 

60  CONTINUE 

CLOSE (12) 

********  subroutine  for  nonlinearity  (discrete  convolution)  ******** 


400 

410 


COMPLEX  FUNCTION  SUM (P,N, NHAR) 
COMPLEX  P (900) 

SUM  =  CMP LX (0 . , 0. ) 

DO  400  M=l,  N-l 

SUM  =  SUM  +  P(M)*P(N-M) 


CONTINUE 

DO  410  M=NHAR,N+1,-1 

SUM  =  SUM  +  2 .  *P  (M)  *C0NJG  (P  (M-N) ) 


CONTINUE 

RETURN 


I 


IV.  SOUND  BEAMS 


SOOMi>  BEAMS 


AJ  a. 


c, 


'jSsUmLt 


_  l*az 

“Z©“  2. 


K?rf<  £<^u<a  Vcrvi  *. 


32-P  rrz* 


x^e.  &  ay- 

2(2  3't’3  2|»  Ce4 


Af-ffracWA  abtorp+icrw  iAo^vUo€a.r((v 

O  Ait  rv:  .or-  e  <->  n  . .  ~  L. J.a  i 


ftssaMp4f<nAS' 


I  BCAr<yec>£  €9aa(rcAi 
re^ovi^c^d  tcty 

-  o 


•  (diCectivfe  lo’eaCfcs) 

•  ^  (K&S  3  (nerfc  ■fcro  dear  .source) 

•  $  £  20°  (net  tea  -far  o+£  ax/s} 

€olaolo6shoj<3  ^fekofoMirO,  ^oVt  flcaiXsi'.  _L£  ( 

3S  U 


320 


GAUSSIAN  BEAMS 


Source  function: 


qo(r )  =  po  exp  [-(r/a)' 


Primary  wave  ( z0  =  ka2/2): 

p0e-a'z 

*l(r’*)  =  w^exp 


>  r 


(r/a)5 


1  — 


Farfield  directivity  (z  >  z0): 

Di(9)  —  exp[—(ka/2)2  tan2  Q\ 


Second  harmonic  pressure: 

j p  e-ct2Z+j(2a\  -at2)zo 

= - T^JTo - exp 


% \<l 


2(r/«): 


1  1-j'z/^oJ 


x  {Ei\j(2ai  —  a2)z0]  -  Ei\j(2ati  -  a2)(z0  -  jz)]} 


Lossless  limit: 


;P0ln(l  -  jz/z0) 


f  \  v - / 

qz(r,z)=  T—^ 


exp 


2(r/a): 


1  -  jz/z0\ 


Farfield: 
q2(0,  z)  = 


— P0lHlfZf2} — tLHD\{6)  exp  {—jkz  tan2 

Zj  Zq 


,  1-0*14.  | -0,1  , 


-inoif 
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lossless  secottb  ge^e^vgkS 

Ri\AdlouM€V\4a(  CciVUpOtA€Vl-t'. 


S«c <yud  (4<arir\fl<nA(c  OsvMftovAe^: 
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UNIFORM  CIRCULAR  SOURCES 


Source  condition: 


qo(r)  =  poH(a-r) 


where  H  is  the  Heaviside  step  function. 


Axial  solution  for  primary  beam  [valid  for  z  ^  a(ka)1^3]: 

qi(Q,z)  =  j2p0s\n(z0/2z)exp(—aiz  —  jz0/2z) 

=  0  z  =  zQ/2mr,  n  =  1,2,... 


*i 


fOT17 


oT  d.oi  0-1 


la 


i/2c 


Farfield  directivity  (valid  for  $  £>  20° ): 


AW  = 


2J\(ka  tan#) 
Ara  tan# 


-is 
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J.  Acoust.  Soc.  Am. 


Averkiou  and  Hamilton 
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cUrer^-'fc.  ^jOd-Vic 
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Figure  4  (Chapter  8) 


Averkiou  and  Hamilton 
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SELF-DEMODULATION  (BERKTAY — 1965) 


Given  the  source  pressure 

P  =  E(t)  sine v0t 

where 


E(t )  =  slowly  varying  modulation  envelope 


Berktay  predicted  that  in  the  farfield  (cr  1)  and  for  strong 
absorption  ( A  >  1),  the  pressure  on  axis  varies  as 


P  oc 


d2E2{t) 
dt 2 


Figure  3  (Chapter  8) 
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FIGURE  12 

PARAMETRIC  TRANSMITTING 
ARRAY  DIRECTIVITY  PATTERNS 
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EXPERIMENT:  fl^«27.3* 
THEORY:  9^0  SO* 
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THEORY:  fl^'21.0* 
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Figure  10-16. -Sum  ftoqoeucy  gutteruR 
retkal  yiIrh  o  —  o,  experimental  wkm 
1  yd,  P2-  85  dire  1  jibar  X  tap*  to 
kHz  (from  Barnard  et  X.  [  16] ). 


for  a 


ittoo  of  48  ft. 


,,f,  *90klte,f,  *101  dire  1  >«fcarX 

k  lecokkf  «i»y,A  *  1*®»  2-8, 5  and  M 


BarearajVxJUlc-fctfjTrucUord)  t  SUcrok-T, 

0-fVS.ft  S2}m3J-  (lUl) 


V.  DISPERSION 


ReLflKflT/O 
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Let 


^  ^  (p\  —  11  i fit? rflal "  CcorJu'^alc 

^  provH'ttes  iMouCrz&opfc  characJesivab'criA.  af  - 


•  cfo-eMie&f  ceacb'cus 

•  Mol-ecular  i/vlora-hVvts 

•  phase  -fraflilfurvis 


ftssu^-e:  #  *-*, 

~33p  s~  n  ~  OU)  (<f  =  relax.at7<r\A  , 

-K  one  ) 


SccffVAdi  order  s4a4€ 


f ' 1  -  d  f 1 +  5B"|V 1  +  'wc‘_  f f£  e 


r^V 


o(i\  c  (£*•') 


OC£2) 


=  dispers(cnA  =  CCt) 


(2cox/m)ai(o)) 
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Frequency  f,  llz 


Figure  IF  13  Log-log  plot  of  sound-absorption  coefficient  versus  frequency  for  sound  in  sir  si 
2<TC  »t  I  atm  pressure  and  with  a  water-vapor  fraction  h  of  4.676  x  10'*  (RH  «  20%).  The  twn 
relaxation  frequencies  are  12,300  Hz  (OJ  and  173  llz  (Nj. 


X  -  X 
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BASIC  EQUATIONS  FOR  BUBBLY  LIQUIDS 


Assume  uniform  distribution  of  identical  bubbles,  each  of  vol¬ 
ume  V  =  Vo  +  v  and  number/ volume  N.  Average  density  p  of 
mixture: 


Pol  p  —  NV  4-  (1  —  NVo)pio/pi 


Expand  densities  (p  =  p0  4-  p',  pi  =  Pio  4-  p\ ,  pg  =  pgo  4-  p'g)  and 
linearize: 


Po 


-Nv 

po4 


Linear  wave  equation: 


(i) 


Rayleigh-Plesset  equation  [V  =  (47t/3)R3]: 

“  3  JLo  j  R  Pg  —  P 

RR  +  -R2  4-  4i/—  =  - 

2  R  Po 

Adiabatic  gas  law: 

Pg/Po  =  (Vo/V^r 

Quadratic  equation  for  bubble  volume: 

i;  4-  6uov  4-  u%v  4-  r/p  =  av2  4-  6(2vi;  4-  i>2)  (2) 

Zaltx>(o'6sfod.<|4  4  Scvl. 


S eccAJt)  H9& Mof^ 

f  =  2  +  CC-. 

/  „  .4CZ\  s  ? 

(y  +■  zi)*?*-9  fate)  f^cyf' 


e4£eef7V)«  coe&ci-e'td  <£  non  I'  (Aeap 


Figure  5.7 
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KdV-BURGERS  EQUATION 


Plane  wave  equation: 


d2p  1  d2p 
dx 2  Cq  dt2 


=  ~PoN 


d2v 
dt 2 


Below  resonance  (a;2  <Wq),  bubble  equation  reduces  to: 


v  = 


1 


6  . 


a 


(jjt 


P - V - «  V  +  ~2V 

1  (jJq  (jJq  LUq 


Volume  terms  are  second  order,  so  use  v  =  —tjp/uJq: 


v  = 


<4 


rf  Srj  rj  .  arj2  2 

-2P+75P+7aP+—P 


Wr 


Substitute  in  wave  equation,  with  Cqq  =  Cq/(1  +  pC): 


d2p  1  d2p 
dx2  Cqq  dt 2 


MW 

UqVq 


d3p  d*p  arj  d2p 2 

W  +  dF  +  Z$HF 


Introduce  slow  scales: 


t'  =  t  —  x/coo  Xi  =  ex 
0(e2)  relation  on  slow  scale: 


dp  ,  dp  u&p  ,d2p 

te=apdt'  +  bdt*  +  cdiz 
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KdV  SOLITONS 


Ignore  losses: 


dp  _  .  dp  ,&p 
Ox  dt'  dt13 


Single  soliton  solution: 


l2  (t  —  x/v\ 
P  —  Po  seen  I — - —  1 


where 


(v  -  Coo)  «  Po' 


Jour  ml  of  Fluid  Mechanic*,  To/.  85,  part  1 


Ficure  13.  Evolution  of  the  inulti*soiiton  pcrturbntion  in  a  liquid  with  He  bubbles. 
<r  =  30;  aflie  =  0*14;  AP,  =  1*74 x  IO»N/m*;  i,  =  0-tSin. 
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m>7? 


Vr  p  :  • 


-T*  I 


;{ 


Fiourr  14.  Evolution  of  a  single  eoliton  in  » liquid  with  He  bubbles, 
i r  a  13*1 ;  trfjl*  a  0*21;  AT,  =  0*028  x  10*Nfm*;  /«  “  0*08  ro. 
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VI.  ACOUSTIC  STREAMING 
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STREAM  fM6-  RARA-METERS 

f  =  f*  +  ft  + 

f  =  f.  +  f ■  +  f  * 

a  =  a,  +  ajc 

OCrt  OCft  <3Cll) 


(jjh-erfi- 


(Ha.66  ptovjJ  Ra.+e: 


(  *<ut  ?/<*£/>  =  f<j>uWs 


=  fc  {  (wac+  ^^)-ds  +  oaV 


define.: 


D  -  3^  ^ 

—  dW<X£!s  -transport  Vi-eloc i4-M 


<-.€  /V!tj  bo  rg (  11  ficcrush'c  Q-br-eam  fl<?  ( 
PbmicaJ  #gcrus,^7cs ,  Acade* ut'd.  Crocs') 
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Co<v\bin>e.  690+7 ii\ud'(<<  ao^  wvoo/i€n+-uw  -fo  eli'M.ivia’te 
<&/\£?+*j.a.n4  +-tm.€  average : 

> i^Uac«^fac-P 

^  p  s=.-^e<^(at*v)ut  +  u((v*Ut)^> 

=  <ArivhV><f  -force  (p-er  ufli+viokme) 


f%  = 


3<tf> 

~ZyT 


ewi^pte* 


Ui  =  a0e  ^ sam («>■+- to*0 

Fx  » otf.aj  e_2#tJt 

=  O  i£  ck—O 


ACOUSTIC  STREAMS ajg- 
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Figure  7 -6. -Acoustic  streaming  for  a  sawtooth  wave 
(from  Statnikov  (13]). 


Figure  7-4.- Acoustic  streaming  from  a  sound  source  in  water.  The  motion  is  made 
visible  by  a  suspension  of  finely  divided  aluminum  (from  Liebermann  ( 1 1  ])• 
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.  1  itflnule. 


Pmc.  Roy.  8 or..  A ,  vol.  KJ4,  PI.  7. 


Ktc.  !h — Ki vi*  sections  of  the  nndo-antinodc  circulation.  Wave-length  o0*o  cm. 


Kio.  UK — Three  sections  of  the  node-nntinode  circulation. 


i 

Kto.  II. — Two  sections  of  the  unde-nutiuodc  circulation. 


I’m*.  12. — t>nc  section,  node  to  nntinodc.  of  the  circulation.  Ware-length  ofl •  o  cm. 


31.  Secondary  streaming  induced  by  an  oscillating  cyl¬ 
inder.  A  long  circular  cylinder  is  oscillated  normal  to  its 
axis  by  a  loudspeaker  in  a  mixture  of  water  and  glycerine. 
Suspended  glass  beads  are  illuminated  in  a  cross  plane  by  a 
stroboscope.  The  amplitude  of  oscillation  is  0.17  of  the 


radius,  and  the  Reynolds  number  based  on  frequency  and 
radius  is  70.  The  steady  second-order  streaming  motion  is 
directed  toward  the  body  along  the  axis  of  oscillation  (indi¬ 
cated  by  arrows)  in  the  inner  region,  ami  opposite  in  the 
outer  region.  Photograph  by  M asakazu  Tatsuno 


VII.  RADIATION  PRESSURE 
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RAi>l#T(6K\ 

[$e&  L<e*U Jo-vM,  2i£j  O^]] 

♦  pH  VV.^t'-f  pj*£id.  •  S"fc'0.d.(  tfVG 


-  £>C  Pressure  P(e\d* 


=  d“t  poLr^cl-e  Melodf^  do<A<?S 


0  fJ 
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Water  deflection  due  to 
D.C.  pressure  at  161  dB. 


Pre-fountain  ripple  at 
161.5  dB. 


Water  fountain  at 
161.5  dB 


Brand.  O.,  and  Freund,  H.  (1934). 
Zeitshrift  Fur  Physik.  92,  385-389. 
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Atmospheric  Refraction  of  Sound 

L  Basics 

Introduction 

Mechanisms 

Wave  equation  (speed  of  sound) 

Reflection 

Ground  models 

Meteorology 

Ray  tracing 

Ray  tracing  lab 


Break 


n.  Advanced 

Problems  with  the  above 
Spherical  wave  reflection 
Fast  field  program 
Residue  series 
Normal  mode 

HL  Graduate 

3-D  FFP,  winds,  density 
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Introduction 

•  Understanding  our  world 

•  Noise 

•  Detection/location 

Research 

•  Solution  oriented 

•  Cyclical 

•  Low  level 


TRANS  2 


Mechanisms  affecting  sound  propagation 

•  Spherical  spreading 

•  Atmospheric  attenuation 

•  Refraction 

•  Diffraction 

•  Ground  absorption 

•  Scattering 

•  Non-linear  effects 
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Atmospheric  sound  propagation 

•  Wavelength  long 

•  Non-linear  effects 

•  Ground  surface  effects 

•  Wind 

•  Meteorological  variability 

•  Turbulence 


TRANS  4 


Wave  Equation 


Continuity 


dp 

at 


+  p0V*v  =  o 


Newton’s  second  law  (Euler  eq) 


Equation  of  state 
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a® 

at 


=  V*<D  using 


Gives 


a2p 

“at2 


For  single  frequencies 

V2p  -  k2p  =  0;  k  =  ® 

Solution  for  k  constant  is  plane  waves: 


P  =  Pe 


i  (k*r  -  tot) 


From  Euler  eq 


TRANS  6 


Reflection  from  ground 


Specific  acoustic  impedance 


Uniform  surface,  plane  waves,  invariance  under 
translation 


=  Zs(co)  =  pc£(co) 
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Ground  Impedance 
Delany  -  Bazeley 
Embleton,  Piercy  and  Daigle 
Attenborough 
Stinson  and  Champoux 
Allard 

Sabatier  and  Bass 

Variables 

Porosity 

Flow  Resistance 
Shape  Factors 
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kj  sin  0j  =  kR  sin  0r 
But  kx  =  kR  =  ®  so  ©i  =0r 
y  dependence  then  -  k  cos  ©j,  k  cos  ©x 
Amplitude  of  reflected  wave  R(@i)  f 
Total  pressure  and  velocity 

p  =  f  e'M  [e-ikyy  +  R(0)i  e^yy 

Vy=  COpC01  f  eik*x  [-er*kyy  +  R(©i)  e^kyy 
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So  at  y  =  o 


1  +  R(0i) 
1  -  R(©i) 


and 


R(0l)  = 


£  (0l)  COS  0J  -  1 
C,  (©i)  COS  ©I  +  1 
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Special  Cases 


Fluid-Fluid  pn,  cn ; 


ci  _  cn 

sin  ©j  sin  ©n 


Pncn 

px  ci  cos(©n) 


Rigid  surface:  Vjn  =  0 
R  =  1 

Pressure  release:  p  =  0 
R  =  -1 


Grazing  Incidence: 
R  =  -1 
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First  Problem: 


Lloyd’s  Mirror 
FM  vs  AM  radio 

Second  problem: 

Is  grazing  same  as  pressure  release? 
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Sound  Speed  not  Uniform 


Temperature  and  relative  humidity  affect  c 

M  =  29(l-h)  +  18h 
y  =  (d+2)/d  =  Cp/cv 

d  =  5(l-h)  +  6h 

Plug  in 

c  =  (331  +  0.6  Tc)(l  +  .16h) 
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Meteorology 

Average  Temperature  Profile: 

Adiabatic  lapse  rate:  -9.8  K/km 
(120  km  Thermosphere) 
Diurnal  Variation 

Heating  Cooling 


Wind  Profiles 

Close  to  logarithmic  near  to  ground 
Decoupled  by  inversions 
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For  the  simple  picture  we  say 

c  =  c(z)  +  w(z)  COS  <j) 


(j)  Angle  between  direction  of  propagation  and 
direction  of  wind  flow 
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RAY  TRACING 


Eikonal  Equation 


Simplest  Case: 


Try 


P  (x,y,z)  =  A(x,y,z)  exp  (ico  [t  -  T  (x,y,z)  /c0]  ) 


A(x,y,z)  Amplitude 


E(x,y,z)  phase 


r(x,y,z)  _  q  is  surface  of  constant  phase 


/ 


*  Vr  in  ray  direction 
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Plug  into  wave  equation 


V2  A(x,y,z)  -  2i©  V  A  •  V  T  /  c0 


-co2  A 


vr  vr 


c0  c0 


2ia>  A^-  +  S^  =  0 

C0  C2 


If  A  and  Vr  slowly  varying  so  that 


I_^<<c^v2r<^ 

A  c2 


CD 

c  c 


n£ 


v2r  <  cq  =  n  co 

c  c  c 


vr<f- n 


TRANS  18 


Then 


tfA^.^A 
C0  Co 


or 


The  Eikonal  Equation 


vr-  vr  =  n2 


n  =  n(x,y,z)  = 


o 


c  (x,y,z) 


Note  Limits  ! 


X 


1  dc 
c  dz 


« 


A  dz 


1 

1 
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Snell’s  Law 


Express  Eikonal  in  terms  of  direction  cosines. 


Express  direction  cosines  in  terms  of  ds. 


Note 


vr 


=  n  so  Eikonal  equation  can  be  written 


a2  +  p2  +  y2  =  1 


where 


(ZEk 

Ivri 


<r 

dx 
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Direction  cosines  in  terms  of  ds 


z 


dz 

ds 


So  combining  solutions  for  7 


dx 

=  n 

1 

dr 

dr 

ds 

n 

dx 

dx 

dy 

=  n 

1 

dr 

dr 

ds 

n 

d\ 

dy 

dz 

=  n 

1 

dr 

dr 

ds 

n 

dz  " 

dz 
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Note: 


dr  dr  n  dr  dr 

=  a  —  +  (3  —  +  y  — 
ds  dx  dy  dz 


so 


=  ana  +  Pnp  +  yny  =  n 
ds 


Taking  ds  of  both  sides  of  each  equation: 


3-D  Snell’s  Laws 


f  M*)  =  - 

ds  V  ds  /  0x 


JL  nJl 
ds  l  ds 


3n 

dy 


-d_fn_dz 

ds  V  ds  > 


3n 

dz 
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n  4^  =  constant  =  n  a0 
ds 

n  4^  =  constant  =  n  B0 
ds 

Choose  x,  z  plane 


z 


dx 


n  cos  0  =  n0  cos  0O 


y  =  Ma  =  sin  0 
ds 
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Radius  of  Curvature 


Put  4^  =  sin©  into  3rd  Snell's  Law 
ds 


-j-  (n  sin  0)  =  4^ 
ds  dz 


use 


dn  _  dn  dz  _ 
ds  dz  ds 


4®-  sin  0 
dz 


d©  _  cos  ©  dn 
ds  n  dz 


cos  Q  dc 
c  dz 


so  curvature 


|d@l  _ 


p  —  c  1 
cos  ©  dc/dz 
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For  linear  sound  speed  change 


R  is  constant  in  layer 

R  -  c  dz 
cos  0  dc 

c  =  q  (1  +  az)  =  cj  +  gz 


p  constant  in  atmosphere,  g  constant  in  layer 
Example:  Huge  -  lm/s/m 

R=  24Q  m/s - L_ —  =  340m 

1  1  m/s/m 
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Construction  of  Rays  for 
Linear  Gradient  Layers 


What  happens? 


Becomes  horizontal  at  z  s.t. 

cos  0  _  1  _  cos  0Q  _ 
c(z)  c(z)  c  (zs) 


c(z)  = 


c(zs) 
cos  (0O) 


If  no  c(z)  > 


c(zs) 

COS  (@o) 


doesn't  come  back 


TRANS  26 


TRANS 


Rays  through  a  Layer: 


©j,  ©f 


calculated  from 


cos  0 
c(z) 


Ar  =  R  |sin@f  -  sin©j 


t  =  R  cos  ©f  -  cos  ©j 


So 

=  jsinQf  -  sin©j| 
|cos0f  -  COS©j 

_  |cOS©f  +  COS@i| 

-  |sin©f  +  sin©; 
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What  about  the  ground? 


© 


We  use  a  mirror  atmosphere. 


TRANS  29 


400 


Area  for  straight  rays  = 

2%  R  cos  ©j  R  d@j 
Area  for  refracted  rays  = 

2k  R  sin  ©f  (4§-)  d©j 

\d0/i 


TRANS  30 


Intensity  is  hteher  if  area  smaller 


focus  factor 


_L  Area  without  refraction 
-L  Area  with  refraction 


2tcRcqs0oR  d0o 

2tiR  (4§-)  d0o  sin© 
\d0/0n 


R  cos0q 
(411)  sin© 

\(1©/q 


What  if 


(dR| 

'd0/o 
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Fermat’s  Principle: 


Rays  travel  in  a  path  which  is  stationary 
with  respect  to  travel  time _ 


Identifies  rays  above  plus 
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PASS  96  Outdoor  Sound  Laboratory 

June  26,  1996 

I.  Draw  rays  coming  from  the  source  using  the  following  rules: 

A  Ray  paths  are  circles  tangent  to  the  initial  direction  -  radii  are  given  on  worksheet 
(+  is  downward  curving;  -  is  upward) 

B.  Angle  of  reflection  =  angle  of  incidence  at  the  ground. 

C.  At  a  layer  interface  the  different  circle  radii  are  tangent  to  the  same  line. 

D.  If  R  =  oo  rays  are  straight  lines. 

II.  Add  additional  ray  paths  as  needed  to  understand  the  ray  field. 

III.  Think  about  what  a  listener  would  hear  at  different  positions. 
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CASE  I 


frg  SN 


CASE  II 


CASE  IV 
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Problems 

CASE  I  Lots  of  energy  close  to  ground 

CASE  II  No  sound  audible  by  ray  trace 

CASE  HI  Lots  of  energy  at  caustics 

CASE  IV  No  sound  heard  since  sound  speed 

never  higher  than  that  at  source 
but _ 

And 


And 


Lloyd’s  Mirror  Effect 


Caustics  -  k/2  phase  change 


TRANS  38 


409 


Computer  Program 

Rays  only  come  back  if  c(z)  >_  c(zs) 

Rays  turn  in  a  layer  if  g(z+)  >  g(z-) 

Ground  treated  as  a  mirror 

and  time  of  flight 

Newton’s  method  used  for  each  class  of  eigenrays 
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Helmholtz  Equation  in  Cylindrical  Coordinates 


Separate  variables  with  an  integral  transform 


p  (K*)  =  (  p(r,z)  J0  (Kr)  rdr 
Jo 

p  (r,z)  =  f  pOM  Jo  (Kr)  KdK 
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Boundary  Conditions 

1 .  Complex  impedance  ground 


2.  Pressure  continuous  at  zs. 
velocity  discontinuous. 

3 .  Pressure  goes  to  zero  at  °° 


Particle 
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Applications 

Spherical  wave  reflection 
Fast  field  program 

Residue  series  solution  for  upward  refraction 

Normal  mode  solution  for  downward 
refraction 
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Spherical  Wave  Reflection 

Point  source,  complex  impedance 
ground 


Uniform  atmosphere  -  e^ 

How  does  this  differ  from  plane  wave 
reflection? 


R(0)  _  cosQ  1 
Zc  COS0  +  1 
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Brekhovskikh 

Converts  to  0  space 
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Solution  by  Steepest  Descents 


Biggest  contribution  at  0O  -  falls  off  rapidly 
Lloyd's  mirror 

x  -  possible  pole  at  Z  cos  0+1=0 
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Spherical  Wave  Reflection  (Donato) 

(Attenborough,  Hayek,  Lawther) 

^  _  e-ikoRi  .  e-ik»R2 

Ri  R2 


2  1 

e-ik0R2 

h 

IN 

OK) 

i 

H* 

OO 

Zccos  ©o  +  1. 

V  rR2 

L  k0R2  J 

+ 


4^  V sin  0n  exp  [iit/4  -  ik0R2  cos  (0P  -  ©0)]. 

jL 7tCl  Z/p 


;  ©p  solution  of  cos  ©p 


JL 

Zc 
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300  Hz,  0.1  m  source,  1 .0  m  receiver 


ui  0‘  V  ©  SP  0‘0  |9Aen 
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range  (m) 
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Break  Up  The  Equation  for  P 


i£  =  -io)poUz 


dU2 

dz 


-i(k2(z)  -  K2)  „  . 

Wo 


2 

copo 


8(z  -  zs 


Two  Ways  To  Solve  For  P(k,z) 


1)  Propagator  Matrix  Method 

2)  Transmission  Line  Analogy 
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Advantages  Of  Each  Approach 


Propagator  Matrix  -  Can  Do  Solid  Layers 

Transmission  Line  -  Can  Have  Reflectionless  Terminations  And 

Numerically  Stable 
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Evaluating  the  integral 

/•  ©O 

p(r,z)  =  I  p  (K,z )  J0  (Kr)  K  d  K 
J  o 

J0  (Kr)  =  1  [Ho1  (Kr)  +  H02(Kr} 

Ho1  (*>  »  VS 

SO 


^  =  V# 


p  (K^r)  eiKr  VK  dK 


This  is  Fourier  transform  of 


P  M  VK 

Use  F.F.T  to  give  p(rm);  m  =  0,  N-l 
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(W)  il|6tbh 
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179  SNVH1 


Wavenumber 
Case  #2 


200 


9S  SNV'dl 


(gp)  UOpCTUSJJV 


Range  (km)  Wind  Speed  (m/s) 

Figure  0.1:  15  Hz  attenuation  excess  (a)  predicted  by  FFP  for  a  simple  model  of  low-level 

wind  shear  (b).  Short  dashes,  upwind.  Long  dashes,  crosswind.  Solid  line,  downwind. 
Atmosphere  assumed  to  be  thermally  homogeneous  ( dT/dz  =  0). 


'64 


l 


O 


400  800  1200 

Rangelml 


1600 


FIG.  6.  Comparison  of  the  data  for  Bondville,  IL  with  coherent  parabolic 
equation  predictions  and  with  the  turbulent  FFP  calculation  f=630  Hr  tat 
a  =  -0.8  m/s,  (b)  .  =  -0.4  m/s,  (c)  .=0  m/s,  (d)  .=0.4  tjs/and Z  0  8 

m/s.  »  are  data  points  with  error  limits;  A— turbulent  FFP  calculation.  The 

continuous  line  is  the  coherent  parabolic  equation  calculation. 
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Questions 


Accuracy  of  FFP  for  upward  refraction 
Role  of  surface  wave  for  upward  refraction 
Criteria  for  downwind  refraction 
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Residue  Series  for  Upward  Propagation 

Point  source,  complex  impedance  ground 
upward  refracting  atmosphere. 

We  pick  a  special  atmosphere 


Equation  for  p  (K,z): 


^E-  +  fki-K2  +  Ml 


dz2 


R 


p  =  o 


Solution  are  Airy  functions 
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BILINEAR  PROFILE 
c  =  c0(l  +^)'1/2 

P  =  -2n  /eiIt/6  AilYx  -  exp^i^-^j 

z>\  [aF(x)  -  q  Ai(x)j  Ai[(x  -  zg)//  exp(i27t/3)] 
1  '  Ai'fx  exp(i27t/3)]  -  q  Ai[x  exp(i2rc/3)] 

l  =  (R/2k02)1'3,  q  =  iko  Z/Zc 
x  =  (k2  -  k02)  P- 

Zeroes  of  Ai'  -  q  Ai  give  Residue  Series. 

q  =  0  zero  of  Ai' 
q  — >  oo  zero  of  Ai 
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Source  m 

\ 

/ — 7 — 7 


General  Behavior 


i 

P  =  W2 


e-ocw  e 


ipw^i>(kw-coO 


P=Im 


-I7t 

*e~6~  bl 


Jeo_\3 

2R2)  J 


TRANS  61 


433 


Residue  Solution:  Upward  Refraction 


TRANS  63 
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Spherical  Wave  <=>  Upward  Refraction  <=>  FFP 

•  Surface  wave  contributes  to  residue  series 

•  Arises  from  zero  of  Airy  function 

•  SWP  Only  important  in  transition  region 

•  FFP  Accurate  to  -120dB 

•  Turbulence  contributes  to  levels  in  the  shadow 
zone 


TRANS  65 


Normal  Mode  Method  for 
Refraction 


Point  source,  complex 
boundary,  downward 
atmosphere. 

Bilinear  profile 


Three  boundary  conditions 


Downward 


impedance 

refracting 
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Normal  Mode  Solution: 


P  (z,k)  =  -2n  ein/6  lAi  (t  +y>)X 


Ai[(x+y<)e  *2jt/3 


Ai'(xei2lt/3)  +  q  Aifrei27[/3) 
Aj(x)  +  q  Ai  tx) 


X  Ai(x+y<) 

/ 


1  =  (R/2ko)1/3 
q  =  ik0l/Zc 
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Residue  Solution:  Downward  Refraction 


swp  -  Subsonic  if  impedance  angle  large 
and  curvature  small. 

7  HI  HcCdVZ  f  fl-M  QaLTS 

'  TRANS  69  1 
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Numerical  Example: 
R  =  200  m 


f=  10  Hz 

Geff  =  300,000  cgs  units 


Impedance  Model: 

Attenborough  low  frequency  approximation 


q  =  .035  +  i  .035 


Re(u) 

-0.984440794E+00 
-0.323741 68 1E+01 
-0.481283475E+01 
-0.615762629E+01 


Im(x) 

0.355343751E-01 
0.1081 13171E-01 
0.727297229E-02 
0.56848055  IE-02 


Re(kn)  Im(kn)  hn 


0.171293223E+00 

0.135401576E+00 

0.103133024E+00 

0.637657608E+01 


0.506772578E-03 

0.195056080E-03 

0.172273623E-03 

0.217787361E-03 


14.0831 

46.3159 

68.8546 

88.0938 
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Magnitude  of  the  spectrum  ,  Mode  strength 


441 


Wavenumber 


ZL  SNVH1 
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Normal  Mode  <=>  Spherical  Wave  <=>  FFP 

•  Impedance  pole  always  contributes 

•  Only  a  true  surface  wave  for  mild  refraction 

•  Normal  mode  strength  and  attenuation  related 
to  the  FFP  spectrum 

•  Low  inversions  don’t  lead  to  large  enhancement 
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Conclusions 

•  Full  wave  needed  at  caustics,  in  shadow  zones, 
and  near  to  ground 

•  Determing  parameter  is  1  =  (R^k,,2)1^3 

*  Surface  wave  behavior  does  not  lead  to 
enhancement  at  long  range 

*  Ideal  cases  help  us  to  understand  FFP  results 
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tbgabr@nadc.navy.mil  State  College,  PA  16804 


Question  #1 


(This  level  is  equivalent  to  an 
1  mg  acceleration  of  0.5  micro-g’s.) 
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Equilibrium  Thermal  Fluctuations: 
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An  atom  in  a  solid  OR  a  mass  on  a  spring  has 
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Sensing  of  Gravity  by  Protozoa 

1.  Fenchel  and  Finlay,  "Geotaxis  in  the  ciliated  protozoon  Loxodes,"  J.  Exp.  Biol.  110,  17-33  (1984) 

2.  Fenchel  and  Finlay,  "The  structure  and  function  of  Muller  vesicles  in  Loxodid  ciliates,"  J.  Protozool.  33,  69-76  (1986) 
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Equilibrium  Thermal  Fluctuations: 
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density  (watts  per  hertz),  or  an  amplitude-squared  per  hertz  (newtons2  per  hertz,  volts2 
per  hertz),  or  an  amplitude  per  root  hertz  (pascals  per  root  hertz,  meters  per  second  per 
root  hertz). 


Noise  Equivalent  Signal: 
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Noise  Equivalent  Sienal: 


Noise  Associated  with  Radiation 
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Noise  Associated  with  Radiation: 
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k-space  density  of  thermal  energy  =  kBT  (L/n) 


Noise  Associated  with  Radiation 
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Equilibrium  Noise  in  a  Geoyhone  (velocity  transducer ); 
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Equivalent  Noise  Generators: 
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Warminster,  pa  18974  Noise-Floor  Measurement 

x4  Array  5/21/95 


Sienal-to-Noise  Ration  A  Useful  Theorem 
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Analysis  of  complicated  circuits  can  often  be  simplified 
by  setting  the  output  load  to  zero  and  calculating  the  ratio 
of  signal  current  to  noise  current. 


S ignal-to -Noise  Ratio  and  Feedback: 


Not  this  way! 


QUIZ :  Question  #2: 
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Shot  Noise: 
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Applies  to  processes  consisting  of  events  that  are: 

(1)  impulse-like 

(2)  independent 


Circuit  for  Basic  Noise  Experiment: 
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1)  Bypass  power  connections 

2)  Use  copper  ground  plane 


Resistor  Noise 

with  and  without  current  flow 


Noise  From  Molecular  Collisions  ( Free-Molecular  Flow): 
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Noise  From  Molecular  Collisions  (Free-Molecular  Flow}: 


Looks  like  a  shot-noise  expression! 
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Molecular-Impact  Noise: 
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Noise  in  Metallic  Conductors 
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Noise  in  Semiconductors: 
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Occupation  (f)  and  Fluctuation 


Occupation  (f)  and  Fluctuation  (df): 
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Noise  as  a  Function  of  Applied  Volta 


Noise  as  a  Function  of  Applied  Voltage: 
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ShotNoise  (Single  Characteristic  Process): 
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Shot  Noise  (Multiple  Characteristic  Processes): 
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Reciprocity: 
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piezoelectric  electron-tunneling 

capacitive  (small-signal) 


Reciprocity  Calibration: 
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It  is  often  inconvenient  and  inaccurate  to  directly  measure 
forces  and  velocities. 


If  a  device  is  reciprocal,  then  the  ratio  of  the  output  potential  to  the 
input  flow  is  the  same  regardless  of  which  port  is  taken  to  be  the  input. 

For  example: 
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Reciprocity  does  NOT  mean  that 
voltage  ratios  are  identical. 
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If  B  does  not  equal  either  C  or  -C,  then  the  device  i 
nonreciprocal. 
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the  fourth  equation  and  allows  solving  for  all  four  responses. 


Recivrocitv  Calibration: 
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Reciprocity  Calibration  ( examples 


Reciprocity  Calibration  (examples 
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Free-Mass  Reciprocity  Calibration 

PCB393A31  SN3001 
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Bessel-Null  Calibration 
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Bessel-Null  Calibration: 
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(where  z,  =  3.83171,7.01559, 10.17347,  13.32369, 


Components  of  a  Sensor  System: 


REVIEW:  HEAT  ENGINES  AND  REFRIGERATORS 
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REVIEW:  SOUND  WAVES 
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Advantages:  reliable,  cheap 

Shortcomings:  big,  inefficient 
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LENGTH  SCALES 

Along  propagation  direction  x: 

Wavelength  A  =  a/f 

Gas  displacement  amplitude  IxJ^IuJ/cj 
Perpendicular  to  x: 

Thermal  penetration  depth  SK  =  \/2K/a;pCp 
Viscous  penetration  depth  6V  =  V2jI/cop 

Relative  sizes: 

=  ME.  <  i 
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’’Audio”  acoustics: 

xj  «  6k«  A 

Thermoacoustic  engines  and  refrigerators: 
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Outline  of  what  will  follow: 
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other  software,  http://rott.esa.lanl.gov/ 
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VISCOSITY  IN  OSCILLATORY  FLOW 

Suppose  acoustic 
oscillation  along  a  wall: 

Momentum  equation: 

^(f"t  +  (V’V)^  =  ~ Vp  +  yUV2V  +  ^+qj^v(v-v) 

+  V/xVv  terms: 

Boundary  condition  v—0  on  surface. 


Linearize  by  using 
P  =  Pm+  p1(x)eitJt 

P  =  Pm(x)  + /01(x,y)ei“t 

V  =  (xUj(x,y)  +  yv1(x,y))eiut 

to  §et  C 

dp  a  % 

lw^ui  =  “  di  +  ^ 

The  solution  is 


U,  »  -  &  ( 1  -  e^(1+i)y/M 

1  icopm  dx  \  x  ^  J 

The  characteristic  dimension  is 

the  viscous  penetration  depth  6U  =  \ /^jx/cop 


DISSIPATION  BY  VISCOSITY  IN  OSCILLATORY  FLOW 
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TEMPERATURE  DUE  TO 
OSCILLATORY  PRESSURE 

Suppose  acoustic 
oscillation  near  a  wall. 

Equation  of  heat  transfer: 
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with  boundary  condition  T=Tm  on  surface. 


Linearize  by  using  p  =  pm  +  Re[pt  eicjt  ] 

T  =  Tm  +Re[T1(x)  e“l] 
P  =  pm  +Re|>1(x)  e™1] 
v  =Re[xu1(x)  eiut] 

d2Tx 

and  get  =  icjpt  +  K^p- 


The  solution  is 


(l+i)x/(5 


The  characteristic  dimension  is  the 

thermal  penetration  depth  6^  =  Y2K/cjpcp 


DISSIPATION  BY  OSCILLATORY  THERMAL  RELAXATION 
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CONNECTION  TO  THINGS  SEEN  ELSEWHERE: 

SOUND  ATTENUATION 
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MOMENTUM  EQUATION 
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ROTT’S  WAVE”  EQUATION 
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How  much  power  is  absorbed  (or  produced) 
in  length  dx  ? 
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ROTT’S  ENERGY  EQUATION 
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SUMMARY,  BASIC  THERMOACOUSTIC  CALCULATION  METHOD 
Numerical  integration  of: 

Momentum  equation  -  gt(  pt ,  <u1>,  Tm) 
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Many  approximate  methods  exist.  Typical  approximations: 

•  fK  and  %  are  given  by  boundary— layer  limit 
°  Tm(x)  is  linear 

•  gas  properties  are  independent  of  temperature 
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CALCULATED  PERFORMANCE 
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5%  conduction  along  stack 
4%  conduction  along  stack 
3%  dT  in  heat  exchangers 


nsider  the  force  exerted  on  a 
here  by  a  moving,  viscous  fluid: 
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the  number  of  independent  units  (m.L.t). 
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Similitude  in  thermoacoustics 


Engine  dimensions:  L,  An{,  hn{,  xjf 
Gas  properties:  y,  a^,  KTe{,  bK,  pK{,  , 

Solid  properties:  Ksj 

SlY)ni\ctY-‘ 

Miscellaneous:  Q,T^Pm, 

Dependent  variables:  /,  T(x,t),  v(x,t),p(x,t),  W. 
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v(x,r)/aref 
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Future:  Details  at  heat  exchanger— stack  interface 

High  amplitude  effects:  turbulence,  streaming.. 
Real  applications 
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About  1/3  of  Carnot  is  typical  of  current  designs. 
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